The Astrophysical Journal, 776:21 (25pp), 2013 October 10 

© 2013. The American Astronomical Society. All rights reserved. Printed in the U.S.A. 


doi : 1 0. 1 08 8/0004-637X/7 76/1/21 


HERSCHEL / PACS SURVEY OF PROTOPLANETARY DISKS IN TAURUS /AURIGA— OBSERVATIONS 
OF [Oi] AND [C ii], AND FAR-INFRARED CONTINUUM 

Christian D. Howard 1,2 , Goran Sandell 1 * * , William D. Vacca 1 , Gaspard Duchene 3,4 , Geoffrey Mathews 5,6 , 
Jean-Charles Augereau 4 , David Barrado 7,8 , William R. F. Dent 9 , Carlos Eiroa 10 , Carol Grady 11,12,13 , Inga Kamp 14 , 
Gwendolyn Meeus 10 , Francois Menard 4,15 , Christophe Pinte 4 , Linda Podio 4 , Pablo Riviere-Marichalar 7 * , 

Aki Roberge 11 * * , Wing-Fai Thi 4 , Silvia Vicente 14 , and Jonathan P. Williams 6 
1 SOFIA-USRA, NASA Ames Research Center, MS 232-12, Building N232, Rm. 146, P.O. Box 1, Moffett Field, CA 94035-0001, USA 
2 Google, 1600 Amphitheatre Parkway, Mountain View, CA 94043, USA 
3 Astronomy Department, University of California, Berkeley, CA 94720-3411, USA 
4 UJF-Grenoble 1/CNRS-INSU, Institut de Planetologie et d’Astrophysique (IPAG) UMR 5274, F-38041 Grenoble, France 
5 Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands 
6 Institute for Astronomy (IfA), University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA 
7 Centro de Astrobiologfa, Depto. Astrofisica (CSIC/INTA), ESAC Campus, P.O. Box 78, E-28691 Villanueva de la Canada, Spain 
8 Calar Alto Observatory, Centro Astronomico Hispano Aleman, C/ Jesus Durban Remon, E-04004 Almerfa, Spain 
9 ALMA SCO, Alonso de Cordova 3107, Vitacura, Santiago, Chile 
10 Dep. de Ffsica Teorica, Fac. de Ciencias, UAM Campus Cantoblanco, E-28049 Madrid, Spain 
11 Exoplanets & Stellar Astrophysics Laboratory, NASA Goddard Space Flight Center, Code 667, Greenbelt, MD 20771, USA 
12 Eureka Scientific, 2452 Delmer, Suite 100, Oakland, CA 96002, USA 
13 Goddard Center for Astrobiology, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
14 Kapteyn Astronomical Institute, Postbus 800, 9700 AV Groningen, The Netherlands 
15 UMI-FCA (UMI 3386: CNRS France, and U de Chile/PUC/U Conception), Santiago, Chile 
Received 2013 March 8; accepted 2013 August 19; published 2013 September 20 


ABSTRACT 

The Herschel Space Observatory was used to observe ~120 pre-main- sequence stars in Taurus as part of the GASPS 
Open Time Key project. Photodetector Array Camera and Spectrometer was used to measure the continuum as 
well as several gas tracers such as [Oi] 63 /tm, [Oi] 145 /rm, [Cn] 158 /rm, OH, H 2 0, and CO. The strongest 
line seen is [O i] at 63 /rm. We find a clear correlation between the strength of the [O i] 63 /rm line and the 63 /rm 
continuum for disk sources. In outflow sources, the line emission can be up to 20 times stronger than in disk sources, 
suggesting that the line emission is dominated by the outflow. The tight correlation seen for disk sources suggests 
that the emission arises from the inner disk (<50 AU) and lower surface layers of the disk where the gas and dust 
are coupled. The [O i] 63 /rm is fainter in transitional stars than in normal Class II disks. Simple spectral energy 
distribution models indicate that the dust responsible for the continuum emission is colder in these disks, leading to 
weaker line emission. [Cn] 158 /rm emission is only detected in strong outflow sources. The observed line ratios 
of [Oi] 63 /rm to [Oi] 145 /rm are in the regime where we are insensitive to the gas-to-dust ratio, neither can we 
discriminate between shock or photodissociation region emission. We detect no Class III object in [O i] 63 /rm and 
only three in continuum, at least one of which is a candidate debris disk. 

Key words: circumstellar matter - planetary systems - protoplanetary disks - stars: pre-main sequence 
Online-only material: color figures 


1. INTRODUCTION 

Low-mass stars are born with protoplanetary disks, canoni- 

cally composed of 1% dust by mass (adopted from the interstel- 
lar medium gas/dust mass ratio), the remaining portion in gas. 
For low-mass stars the gas has largely disappeared on timescales 

of 6-7 Myr (Haisch et al. 2001; Hernandez et al. 2008). Under- 

standing the early evolution of stars requires the understanding 

of their accompanying protoplanetary disks of gas and dust. 

Gas in Protoplanetary Systems (GASPS) is a large Open Time 

Key Project on the Herschel Space Observatory (Pilbratt et al. 

2010) studying the evolution of gas in protoplanetary disks. 

In total the project surveyed ~250 nearby (^200 pc) low- and 
intermediate-mass stars (0.3-8 M©), from young (0.5 Myr) stars 

with massive disks to older (30 Myr) stars with very little dust 

excess. For this study we are using the Photodetector Array 

Camera and Spectrometer (PACS; Poglitsch et al. 2010) to 
observe the fine-structure lines of [Cn] at 158 /rm and [Oi] 

at 63.2 and 145 /rm as well as several H 2 0, OH, and high 
rotational transitions of CO tracing mainly hot gas. We also 


use the PACS imager to measure broadband continuum fluxes 
for a large portion of the sample. The goal is to understand the 
transition from gas rich protoplanetary disks to gas poor debris 
disks, and the timescales on which the gas dissipates. For a more 
detailed description of GASPS, see Dent et al. (2013). 

The largest sample of GASPS sources are in the Taurus star 
forming region. As one of the closest star forming regions to the 
Sun (distance of ~ 140 pc), the Taurus group provides a unique 
opportunity to obtain a statistical sample of largely coeval 
protoplanetary disks encompassing a wide range of disk masses. 
The majority of GASPS sources in Taurus are class II objects 
(classical T Tauri stars, CTTSs, with spectral types ranging 
from mid-F to early M, and with stellar masses of ~0.3 to 
1.4 M q ), having dissipated their envelope of accreting material 
and leaving behind the protostar and a disk of gas and dust (see, 
e.g., Kenyon et al. 2008). Of all the gas probes used by GASPS, 
the [O i] 63 jim line is by far the most sensitive gas tracer, at 
least for protoplanetary disks (Pinte et al. 2010). 

The [O i] 63 /xm line is known to be strong in photodisso- 
ciation regions (PDR; Tielens & Hollenbach 1985), and also 
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in shocks (Hollenbach & McKee 1989; Flower & Pineaut des 
Forets 2010). Because of this, it is necessary to isolate the com- 
ponent of the emission arising from the disk itself from that 
arising from jets or outflows, which are common in young T 
Tauri stars. The velocity resolution of PACS is insufficient to 
separate the low velocity disk emission from the outflow emis- 
sion and for most objects the outflow is spatially unresolved 
in [Oi] 63 /xm. Podio et al. (2012) have spatially resolved the 
[O i] 63 /xm emission in a few of the strongest jet sources in the 
Taurus sample. In our analysis we therefore distinguish between 
sources that are already known to have jets/outflows based on 
previous observations and sources which have very weak or un- 
detectable outflows; in the latter the [O i] emission, if detected, is 
expected to be dominated by emission from the protoplanetary 
disk. The relative importance of FUV and X-ray irradiation for 
the [O i] line emission from the Taurus disk sources is explored 
in a forthcoming paper by Aresu et al. (2013). Detailed SED 
modeling and molecular emission lines will be discussed in a 
future paper. 

In Section 2 we discuss the Taurus sample, In Section 3 we 
discuss the observations and reduction for both spectroscopy 
and photometry as well as pointing issues that arose in some 
observations. Section 4 discusses the results for the [O i] and 
[C n] fine structure lines, followed with a discussion of what 
we learned in Section 5. Section 5.6 gives a brief summary of 
our results for Class III objects and in Section 5.7 we comment 
on a few individual objects. Section 6 gives the summary and 
conclusions of this study. 

2. THE TAURUS-AURIGA DARK YOUNG 
STELLAR ASSOCIATION 

The Taurus-Auriga cloud complex is one of the closest 
(140 pc; Bertout et al. 1999; Torres et al. 2009) active star 
forming regions, with a young stellar association containing at 
least 250 cluster members (Rebull et al. 2010, 2011). It does 
not have any O and B stars, but it has a rich population of 
pre-main- sequence (PMS) low-mass T Tauri stars and brown 
dwarfs. The Taurus-Auriga dark cloud complexes extend over 
~10° of the sky with a depth of ~20 pc (Torres et al. 2009). 
The youngest stars, Class I protostars, are found toward the 
center of opaque dark clouds. Likewise the CTTSs (Class II) 
are also concentrated toward dark clouds, whereas the Class III 
objects, weak-lined T Tauri stars, are more widely distributed 
over the whole region. Our Taurus sample also includes nine 
confirmed or suspected transitional disk objects. Transitional 
disks, which were first discovered by Strom et al. (1989), are 
believed to be transition objects between CTTSs and weak-lined 
T Tauri stars. They have large mid- to far-infrared excesses, but 
weak (pre-transitional disks) or no excess in the near-infrared 
(Najita et al. 2007; Espaillat et al. 2011). Both pre-transitional 
and transitional disks have large gaps or cavities largely void 
of gas and dust in the inner disk, which in several cases have 
been resolved with millimeter interferometry (Pietu et al. 2006; 
Dutrey et al. 2008; Brown et al. 2009; Andrews et al. 2011). 
The majority of the PMS population in Taurus has an age of 
1-3 Myr, although there is certainly an age spread. Some stars 
are as young as 0.5 Myr and there are a few objects as old as 
15-20 Myr (Hartmann 2003; Kiigiik & Akkaya 2010). 

Some of the younger, more deeply embedded T Tauri stars 
drive Herbig-Haro (HH) flows or well-collimated jets, typi- 
cally seen in the optical [Oi], Ha, [Nn] AA6548, 6583, and 
[S n] AA6717, 6731 emission lines (Mundt & Fried 1983; Gomez 
et al. 1997). Since we try to discriminate between disk emission 


and outflow emission in our GASPS observations, it is there- 
fore important to identify which stars power optical jets or HH 
flows. For this we have mainly used the tabulation provided by 
Kenyon et al. (2008), but we have also identified a few additional 
ones: AA Tau, DL Tau, and V773 Tau, see Section 5.7, IRAS 
04385+2805 (Glauser et al. 2008), and SU Aur (Chakraborty 
& Ge 2004). Stars that we define as jet/outflow sources must 
have a jet imaged in Ha, [O i] A6300, [S n] A6371 or be associ- 
ated with HH objects, show a high velocity molecular outflow, 
or a broad (>50 km s -1 ), typically blue-shifted, emission line 
profiles in [Oi] A.6300 (see, e.g., Hartigan et al. 1995). 

For the most part, our analysis leaves aside the issue of 
multiplicity. About half of our targets are known binary or 
higher-order multiple systems, with separation ranging from 
<0.1 AU to as much as 2000 AU. On an individual object 
basis, the interpretation of our observations can be severely 
affected by the possible presence of several circumstellar and/or 
circumbinary disks. However, as discussed in Section 5.5, it 
appears that the effect of multiplicity on both far-infrared 
continuum and line emission is relatively modest. 

3. OBSERVATIONS AND DATA REDUCTION 

The observations reported herein were carried out with 
Herschel from 2010 February through 2011 March. We obtained 
photometry and spectroscopy utilizing the PACS instrument 
aboard Herschel. For a description of PACS observing modes, 
see Poglitsch et al. (2010). 

3.1. Photometry 

All photometric observations of our targets in Taurus were 
made in “mini-scan” mode, i.e., by doing short, 3' scans with 
medium scan speed of 20" s -1 , and with 10 4" steps in the 
direction orthogonal to the scan. Executing such an observation 
takes 276 s. Two bands, either blue and red, or green and 
red, are recorded simultaneously. The blue filter is centered 
on 70 /xm, the green filter is centered at 100 /xm and the red at 
160 /xm (see, e.g., Poglitsch et al. 2010). All targets have been 
observed in the 70 /xm and the 160 /xm bands and many were 
also observed in the 100 /xm filter (see Table 1). Almost all the 
Taurus photometric observations were carried out in two scan 
directions (70° and 110°) providing a smoother background and 
hence more accurate photometry. However, some of the early 
photometry was done with only a single blue/red mini-scan. 
Only T Tauri is spatially resolved at 70 /xm. 16 An extended disk 
is, however, detected around HD 283759, an F5 star, but proper 
motion studies (Massarotti et al. 2005) show conclusively that 
the star is not a member of the Taurus association. 

All the photometric data were originally reduced using the 
Herschel Interactive Processing Environment, HIPE 17 v7.0, 
but we later re-reduced all detected sources in HIPE v9.0, 
where the reduction scripts provide iterative masked high 
pass filtering, resulting in better deglitching and smoother 
backgrounds, especially in the red. For the blue filter the results 
were essentially identical, i.e., within l%-2%. The final images 
were imported to the STARLINK program Gaia for photometry 
and further analysis. For strong sources we performed aperture 


16 The emission from UZ Tau is also extended, but here the emission 
originates from two binaries UZ Tau E and UZ Tau W, which are separated by 
376 and which both are surrounded by disks, see, e.g., Harris et al. (2012). 

17 HIPE is a joint development by the Herschel Science Ground Segment 
Consortium, consisting of ESA, the NASA Herschel Science Center, and the 
HIFI, PACS, and SPIRE consortia. 
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Table 1 

GASPS Taurus Observations: PACS Photometry 


Name 

SED 

Class 

70 jLim 
(Jy) 

100 /urn 

(Jy) 

160 /urn 

(Jy) 

Comments 

CIDA2 

III 

<0.009 

<0.009 

<0.022 

Empty field 

CIDA 10 

III 

<0.008 


<0.028 

Empty field 

CoKu Tau 2 

I 

1.30 ±0.20 

1.32 ±0.20 

1.53 ±0.30 

34" west of HL Tau 

CoKu Tau 4 a 

II 

1.15 ±0.06 

0.85 ± 0.08 

0.64 ±0.10 


CW Tau b 

II 

1.79 ±0.18 


1.82 ±0.18 

Same field as V773 Tau 

CX Tau c 

T 

0.33 ± 0.03 

0.25 ± 0.03 

0.25 ± 0.03 

SI 

CYTau 

II 

0.22 ± 0.01 

0.29 ± 0.03 

0.35 ± 0.04 

SI 

DG Tau 

II 

23.8 ± 2.4 

23.3 ± 2.3 

18.30 ±2.0 

Cloud emission, S2 

DG Tau B 

I 

15.4 ± 1.5 

18.0 ± 1.8 

18.10 ±2.0 

Same field as DG Tau 

DH Tau A 

II 

0.47 ± 0.05 

0.48 ± 0.05 

0.41 ± 0.04 


DI Tau AB 

III 

<0.01 

<0.01 

<0.03 

Near DH Tau 

DK Tau A 

II 

1.17 ±0.12 

1.07 ±0.11 

0.83 ± 0.08 


DL Tau 

II 

1.32 ±0.13 

1.54 ±0.15 

1.93 ±0.19 

SI 

DM Tau 

T 

0.78 ± 0.08 

0.85 ± 0.08 

0.78 ± 0.08 


DO Tau d 

II 

5.1± 0.80 

4.7± 0.70 

4.1± 0.70 


DP Tau 

II 

0.67 ± 0.07 

0.36 ± 0.04 

0.20 ± 0.03 

Cloud emission 

DQ Tau e 

II 

1.37 ±0.10 

1.28 ±0.12 

1.02 ±0.10 

SI, field includes Haro 6-37 

DS Tau 

II 

0.22 ± 0.02 

0.33 ± 0.03 

0.26 ± 0.03 


FF Tau AB 

III 

<0.009 

<0.010 

<0.028 

S2 

FM Tau b 

II 

0.47 ± 0.07 


0.24 ± 0.03 


FO Tau AB 

T 

0.51 ±0.05 

0.39 ± 0.04 

0.21 ± 0.03 

SI 

FP Tau c 

II 

0.33 ± 0.03 

0.34 ± 0.06 

0.38 ± 0.04 

SI 

FT Tau 

II 

0.73 ± 0.07 

0.95 ± 0.09 

1.27 ±0.19 


FW Tau ABC 

III 

0.030 ± 0.004 

0.033 ± 0.004 

0.070 ± 0.040 

SI, cloud emission 

FX Tau A 

II 

0.39 ± 0.04 

0.36 ± 0.04 

0.22 ± 0.03 


GG Tau Aab 

II 

4.22 ± 0.42 

6.56 ± 0.66 

8.41 ± 0.84 

Extended 

GG Tau Bab 

II 

0.21 ± 0.02 




GH Tau AB 

II 

0.37 ± 0.04 

0.52 ± 0.05 

0.27 ± 0.03 

SI, field includes V807 Tau 

GI Tau f 

II 

0.67 ± 0.07 

0.50 ± 0.05 

0.33 ± 0.06 


GK Tau f 

II 

0.90 ± 0.09 

0.59 ± 0.06 

0.35 ± 0.06 


Haro 6-37 B e 

II 

0.96 ±0.10 

0.90 ± 0.09 

1.25 ±0.12 


HBC 347 

III 

<0.009 

<0.009 

<0.019 

S2 

HBC 352/353 

III 

<0.007 


<0.017 

SI 

HBC 354/355 

III 

<0.007 


<0.022 

Empty field 

HBC 356/357 

III 

<0.010 

<0.011 

<0.023 

SI 

HBC 358/359 

III 

<0.007 


<0.020 

S2 

HBC 360/361 

III 

<0.009 

<0.010 

<0.020 

Empty field 

HBC 362 

III 

<0.008 


<0.019 

S4 

HBC 372 

III 

<0.009 

<0.009 

<0.021 

Empty field 

HBC 376 

III 

<0.009 

<0.009 

<0.016 

S3 

HBC 388 

III 

<0.007 


<0.018 

Empty field 

HBC 392 

III 

<0.007 


<0.013 

Empty field 

HBC 407 

III 

<0.009 

<0.009 

<0.017 

Empty field 

HBC 412 AB 

III 

<0.009 

<0.009 

<0.018 

Empty field 

HD 283572 

III 

<0.008 

<0.009 

<0.027 

SI 

HD 283759 s 

S 

0.52 ± 0.05 

0.69 ± 0.07 

0.57 ± 0.06 

SI 

HL Tau 

I 

75.3 ± 7.5 

72.3 ± 7.2 

54.7 ± 8.2 

Cloud emission, 4 T Tauri stars, SI 

HN Tau A 

II 

1.02 ±0.10 


0.61 ± 0.06 

S4 

HO Tau 

II 

0.10 ±0.01 

0.11 ±0.01 

0.14 ±0.02 

S2 

HV Tau C d 

I? 

1.55 ±0.16 

1.68 ±0.18 

1.45 ±0.24 


IQ Tau 

II 

0.74 ± 0.07 

0.78 ± 0.08 

0.81 ± 0.08 

S2 

IW Tau 

III 

<0.009 

<0.010 

<0.030 

SI 

J043 05 1 7 1 +244 1475 

II 

2.85 ± 0.29 

3.11 ±0.31 

3.06 ±0.31 

Cloud emission, 35" north of ZZ Ta 

Jl-4872 

III 

<0.009 

<0.009 

<0.027 

S3 

Jl-507 

III 

<0.008 

<0.008 

<0.021 

Empty field 

Jl-665 

III 

<0.009 

<0.009 

<0.025 

Empty field 

J2-2041 

III 

<0.008 

<0.009 

<0.013 

SI 

JH 108 

III 

<0.008 

<0.009 

<0.021 

S2 

JH 223 

II 

0.114 ±0.011 

0.109 ±0.011 

0.076 ± 0.016 


JH 56 

II 

0.027 ± 0.003 

0.017 ± 0.003 

<0.04 

Cloud emission 

L 1551-51 

III 

<0.009 

<0.009 

<0.016 

Empty field 

L 1551-55 

III 

<0.009 

<0.009 

<0.016 

Empty field 

LkCal 

III 

<0.009 

<0.009 

<0.021 

Empty field 

LkCa 3 

III 

<0.009 

<0.009 

<0.017 

SI 

LkCa 4 

III 

<0.009 

<0.009 

<0.034 

Empty field 
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Name 

SED 

Class 

70 /xm 
(Jy) 

100 /xm 

(Jy) 

160 /xm 

(Jy) 

Comments 

LkCa 5 

III 

<0.009 

<0.009 

<0.029 

SI, cloud emission 

LkCa 14 

III 

<0.009 

<0.009 

<0.028 

S2 

LkCa 15 

T 

1.23 ±0.12 

1.48 ±0.15 

1.78 ±0.18 


LkCa 19 

III 

<0.009 

<0.009 

<0.015 


LkCa 21 

III 

<0.012 


<0.051 

In same field as RY Tau 

RW Aur A 

II 

2.47 ±0.15 

2.94 ± 0.29 

1.74 ±0.17 

S2 

RYTau 

T? 

14.13 ± 1.40 


8.81 ±0.88 

Cloud emission 

SAO 76411 

Ill 

<0.012 

<0.011 

<0.020 mJy SI 


SAO 76428 

III 

<0.009 

<0.009 

<0.020 mJy S2 


StHA 34 h 

II 

0.053 ± 0.006 

0.039 ± 0.005 

0.025 ± 0.025 

S2 

Tau L1495 14 

I 

1.11 ±0.22 


1.45 ±0.29 

In same field as V773 Tau, edge of field 

T Tau NS 1 

n/i 

170 ± 20 

143 ± 21 

100 ± 20 

Extended, cloud emission 

UZ Tau Eab+Wab 

ii 

2.03 ± 0.20 

2.13 ±0.21 

2.04 ± 0.20 


V397 Aur 

hi 

0.016 ± 0.003 


0.016 ± 0.007 


V773 Tau b 

ii 

0.81 ±0.08 


0.25 ± 0.03 


V807 Tau AB 

ii 

0.51 ±0.05 

0.72 ± 0.07 

0.35 ± 0.04 

S2, field includes GH Tau 

V819 Tau 

ii/iii 

0.030 ± 0.003 

0.017 ± 0.003 

<0.050 

SI, faint cloud emission 

V826 Tau AB 

hi 

<0.009 

<0.009 

<0.017 

Empty field 

V827 Tau 

hi 

<0.009 

<0.009 

<0.018 

Empty field 

V830 Tau 

hi 

<0.009 

<0.010 

<0.030 

Cloud emission 

V836 Tau 

T 

0.21 ± 0.02 

0.30 ± 0.03 

0.27 ± 0.05 

SI 

V927 Tau 

III 

<0.009 

<0.009 

<0.032 

Empty field 

V928 Tau 

III 

<0.009 

<0.011 

<0.053 

Cloud emission 

VI 096 Tau 

III 

<0.009 

<0.010 

<0.030 

Cloud emission. 

V1213 Tau 

I 

0.41 ± 0.04 

0.51 ±0.05 

0.59 ± 0.012 

In same field as HL Tau 

VY Tau AB 

II 

0.30 ± 0.03 

0.29 ± 0.03 

0.22 ± 0.02 

S2 

Wa Taul 

III 

<0.009 

<0.009 

<0.014 

S2 

XZ Tau AB 

II 

5.77 ± 0.87 

4.71 ± 0.71 

4.13 ±0.80 

23" east of HL Tau 

ZZ Tau 

II 

0.052 ± 0.005 

0.051 ±0.005 

0.025 ± 0.02 



Notes. Upper limits are 3cr. The commonly used source notations CIDA NN and Jl-NN are not know to SIMBAD, which instead uses [BCG93] NN 
and [HJS91] NN, respectively. In the comment field Sn ( n = 1, 2, 3 etc.) lists the number of “background” sources we find in each field. 
a Star surrounded by extended nebulosity, see text. 
b V773 Tau, FM Tau, CW Tau and Tau LI 495 14 in the same field. 
c CX Tau and FP Tau in the same field. 

d DO Tau and HV Tau in the same field. Cloud emission, see text. 
e Haro 6-37 and DQ Tau in the same field. 

f GI Tau and GK Tau in the same field; some nebulosity to the east in red. 
g Emission is extended! Not a member of the Taurus association (Massarotti et al. 2005). 
h Stronger source within 1 8" to the NW. 

1 Extended. 


photometry with a 12" aperture, while we used 6"-9" apertures 
for faint sources and applied the aperture corrections determined 
by the PACS team (Muller et al. 2011). An annulus extending 
from 60"-70" was used for estimating the sky background. In 
some cases, where stars were near the edge of the image, where 
we had multiple stars, or where there was extended background 
emission in the image, the sky emission was estimated from the 
average of several, “clean” areas in the image. For the 70 /xm 
and 100 /xm bands the background is very close to zero. In a few 
cases the star was surrounded by emission from the surrounding 
cloud. This emission is mostly seen at 160 /xm, but can be 
visible in the 100 /xm band. In these cases we used a 6"-8" 
aperture, with the appropriate aperture correction applied. The 
point source calibration for the PACS images has an absolute 
accuracy of 3% in the 70 /xm and 100 /xm band and better than 
5% at 160 /xm (Muller et al. 2011). In our analysis we do not 
reach such accuracy. Many of our target stars are embedded 
in or are in the direction of dark clouds, and variations in 
the background limits the photometric accuracy. We therefore 


estimate our photometric accuracy to be about 10%. In the 
70 /xm and 100 /xm bands the 3 a upper limits are typically about 
9 mJy. In the 160 /xm band the noise is much higher due to the 
presence of “cirrus”-like emission from dust clouds. Typical 3a 
upper limits are about 20 mJy. All flux densities are given in 
Table 1 as measured, i.e., without color corrections. The color 
corrections, however, are small (<2%) for all our targets and do 
not affect our analysis. We also list the SED classes in Tables 1 
and 2. The SED classification comes primarily from Luhman 
et al. (2010) and or Rebull et al. (2010); if they differ we list 
both. Some weak-line T Tauri star (WTTS) were not observed 
by Luhman et al. or Rebull et al. If we found no indication for 
an infrared excess we assigned them as Class III. Transitional 
disks are labeled T. 

Some fields also show faint unrelated sources. Almost all of 
them are extragalactic background sources. A few, however, do 
coincide with faint 2MASS and WISE sources, and could be 
low-mass T Tauri stars. We therefore also list the number of 
“background” sources we find in each field in Table 1. 
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Table 2 

PACS [Oi] 63, [Oi] 145, and [Cn] 158 Flux Densities and Other Possible Detections 


Name 

Jet/ 

Outflow 

SED 

Class 

[Oi] 63 

(10“ 16 Wm“ 2 ) 

[Oi] 145 
(10“ 16 Wm“ 2 ) 

[Cn] 158 
(10“ 16 Wm“ 2 ) 

Other Fines 

AA Tau 

Y 

II 

0.22 ± 0.02 

0.02 ± 0.01 (<0.02) 

<0.02 

1,2, 4,5,9 

BP Tau 

N 

II 

0.10 ±0.03 



1 

Cl Tau 

N 

II 

0.33 ± 0.05 

0.02 ± 0.02 (<0.07) 

<0.07 

9 

CIDA2 

N 

III 

<0.08 




CoKu Tau/4 

N 

II 

0.22 ± 0.02 




CW Tau 

Y 

II 

0.72 ± 0.04 

0.05 ± 0.02 

0.05 ± 0.01 (<0.05) 

9 

CX Tau 

N 

T 

0.07 ± 0.03 (<0.08) 




CYTau 

N 

II 

0.12 ±0.04 




DE Tau 

N 

II 

0.07 ± 0.06 (<0.12) 

0.04 ± 0.02 

<0.04 

9 

DF Tau 

Y 

II 

0.61 ± 0.06 

0.06 ± 0.04 (<0.09) 

0.05 ± 0.02 (<0.06) 


DG Tau 

Y 

II 

13.40 ±0.17 

0.84 ± 0.05 

1.87 ±0.04 

4, 5, 7, 8, 9 

DG Tau B 

Y 

I 

4.26 ±0.12 

0.27 ± 0.03 

0.13 ±0.02 

9 

DHTau 

N 

II 

0.07 ± 0.02 (<0.10) 




DITau 

N 

III 

<0.14 




DK Tau 

N 

II 

0.16 ±0.03 

<0.10 

<0.15 


DL Tau 

Y 

II 

0.22 ± 0.02 

<0.14 

<0.09 

1 

DM Tau 

N 

T 

0.07 ± 0.02 

<0.06 

<0.06 

1 

DN Tau 

N 

T 

0.06 ± 0.02 

<0.06 

0.05 ± 0.02 (<0.09) 


DO Tau 

Y 

II 

0.71 ±0.10 

0.08 ± 0.01 

0.03 ± 0.01 (<0.03) 

2,9 

DP Tau 

Y 

II 

1.48 ±0.13 

0.07 ± 0.02 

0.12 ±0.02 

5,9 

DQ Tau 

N 

II 

0.21 ± 0.04 

<0.04 

0.03 ± 0.01 (<0.04) 


DS Tau 

N 

II 

0.09 ± 0.02 




FF Tau 

N 

III 

<0.10 




FM Tau 

N 

II 

0.10 ± 0.02 (<0.13) 

<0.06 

<0.06 


FO Tau 

N 

T 

0.12 ± 0.05 (<0.12) 

<0.05 

<0.04 


FQ Tau 

N 

II 

<0.09 

<0.03 

<0.06 


FS Tau A 

Y 

II /Flat 

3.58 ± 0.05 

0.18 ±0.01 

0.16 ±0.01 

1,4, 5, 7, 8, 9, 10, 11 

FT Tau 

N 

II 

0.17 ±0.05 




FW Tau 

N 

III 

0.04 ± 0.02 (<0.08) 




FX Tau 

N 

II 

<0.14 




GG Tau Aab 

N 

II 

0.51 ±0.04 

0.03 ± 0.01 (<0.05) 

<0.04 

4, 5,9 

GH Tau 

N 

II 

0.08 ± 0.04 (<0.16) 

<0.03 

<0.02 


GI/GK Tau a 

N/Y 

II 

0.31 ±0.14 




GM Aur 

N 

T 

0.24 ± 0.05 

<0.06 

<0.10 


GO Tau 

N 

II 

<0.06 




Haro 6-5B 

Y 

I 

0.68 ±0.11 

0.04 ± 0.01 (<0.04) 

0.06 ± 0.01 


Haro 6-13 

Y 

II 

0.70 ± 0.05 

0.03 ± 0.01 (<0.04) 

<0.05 

2, 4, 5, 9 

Haro 6-37 

N 

II 

0.10 ±0.02 




HBC 347 

N 

III 

<0.12 




HBC 356 

N 

III 

<0.08 




HBC 358 

N 

III 

<0.14 




HD 283572 

N 

III 

<0.08 




HK Tau 

N 

II 

0.34 ± 0.02 

<0.04 

<0.04 

9 

HF Tau 

Y 

I 

5.13 ±0.05 

0.28 ± 0.02 

0.36 ± 0.04 

1,2, 3, 4, 5, 7, 8,9 

HN Tau A 

Y 

II 

0.41 ± 0.02 

0.04 ± 0.01 

<0.04 

1 

HO Tau 

N 

II 

<0.10 




HV Tau C 

Y 

I? 

0.52 ± 0.08 




IP Tau 

N 

T 

0.06 ± 0.02 (<0.07) 




IQ Tau 

N 

II 

0.15 ±0.03 

<0.11 

<0.07 


IRAS 04158+2805 

Y 

I 

0.57 ± 0.04 

0.05 ± 0.02 

<0.08 


IRAS 04385+2550 

Y 

II 

0.49 ± 0.02 




Jl-4827 

N 

III 

<0.09 




FkCal 

N 

III 

<0.08 




FkCa 3 

N 

III 

<0.10 




FkCa 4 

N 

III 

<0.10 




FkCa 5 

N 

III 

<0.10 




FkCa 7 

N 

III 

<0.09 




FkCa 15 

N 

T 

0.10 ±0.02 

<0.09 

<0.11 


RW Aur 

Y 

II 

1.54 ±0.05 

0.12 ±0.03 

<0.10 

9, 10 

RYTau 

Y 

T? 

1.05 ±0.05 

0.07 ± 0.02 

<0.10 

1,2,9 

SU Aur 

Y 

II 

0.86 ± 0.03 

0.06 ± 0.02 (<0.06) 

0.07 ± 0.02 

9 

T Tau 

Y 

H/I 

81.80 ±0.31 

4.06 ± 0.04 

2.29 ± 0.02 

1,2, 3, 4, 5, 6, 7, 8,9, 10, 11 

UX Tau A 

N 

T 

0.34 ± 0.02 




UY Aur 

Y 

II 

3.14 ±0.04 

0.17 ± 0.01 

0.16 ±0.01 

1,2, 4,5,7, 8, 9, 10 
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Name 

Jet / 

Outflow 

SED 

Class 

[Oi] 63 

(10 — 16 W m -2 ) 

[Oi] 145 
(10 -16 W m -2 ) 

[Cn] 158 
(10 “ 16 Wm" 2 ) 

Other Lines 

UZTau 

Y 

II 

0.45 ±0.14 

< 0.11 

<0.07 


V710 Tau 

N 

II 

0.10 ±0.06 (< 0 . 12 ) 




V773 Tau 

Y 

II 

0.65 ± 0.03 



1 

V807 Tau 

N 

II 

0.14 ±0.05 

<0.05 

0.02 ± 0.01 (<0.04) 


V819 Tau 

N 

II/III 

<0.09 




V836 Tau 

N 

T 

0.05 ± 0.04 (<0.06) 




V927 Tau 

N 

III 

<0.08 




VI 096 Tau 

N 

III 

< 0.12 




VY Tau 

N 

II 

< 0.11 




XZ Tau 

Y 

II 

3.61 ± 0.09 

0.18 ±0.04 

0.24 ± 0.02 

1,2, 4,5,7, 8,9 

ZZ Tau 

N 

II 

< 0.10 





Notes. For sources where no line is detected, the 3cr upper limit is reported. Possible detections of other emission lines determined through visual 
inspection. The numbers refer to 

(1) O-H2O 818 -»• 7 07 

(2) CH + 7 = 5 -» 4 

(3) CO 7 = 36 -> 35 

(4) 0 -H 2 O 423 3 12 

(5) OH l/2-3/2hfs 

(6) CO 7 = 33 -> 32 

(7) p-H 2 0 3 2 2 ^2 11 

(8) CO 7 = 29 28 

(9) CO 7 = 18 -> 17 

(10) 0 -H 2 O 212 -> lot 

(11) o-H 2 0 2 2 i -> 2i2 


a Gi Tau and GK Tau are unresolved in spectroscopy. 

Table 3 

PACS Spectroscopy Settings 

Obs. Mode Setting 

Grating Order 

Camera 

Observed Range 
(/xm) 

Species 

Transition 

Wavelength 

(/xm) 

LineSpec A 

3 

Blue 

62.68-63.68 

[O 1] 

3 Pi -^ 3 P 2 

63.184 





0 -H 2 O 

818 7 07 

63.324 


1 

Red 

188.77-190.30 

DCO + 

7 = 22 -> 21 

189.570 

RangeSpec B 

2 

Blue 

72.00-73.05 

CH + 

7 = 5^4 

72.14 





CO 

7 = 36 -> 35 

72.843 


1 

Red 

144-146.1 

CO 

7 = 18 -> 17 

144.784 





[O 1] 

3 Po Pi 

145.525 

C 

2 

Blue 

78.55-79.45 

0 -H 2 O 

423 — ^ 3i2 

78.741 





OH 

5 — > § hfs 

79.11/79.18 





CO 

7 = 33 -> 32 

79.360 


1 

Red 

157.1-158.9 

[C 11] 

2 P 3 /2 < Pl/2 

157.741 





p-H 2 0 

331 4 0 4 

158.309 

D 

2 

Blue 

89.45-90.50 

p-H 2 0 

322 2 n 

89.988 





CH+ 

7 = 4^3 

90.02 





CO 

7 = 29 -> 28 

90.163 


1 

Red 

178.9-181.0 

0 -H 2 O 

2 l 2 lot 

179.527 





CH+ 

J = 2^1 

179.610 





0 -H 2 O 

221 212 

180.488 


Note. The primary lines are highlighted in bold. 


3.2. Spectroscopy 

3.2.1. Observations 

We observed the Taurus sample with PACS using both 
chopped line and range spectroscopy modes, targeting primarily 
[Oi] 63 and 145 /xm and the [Cn] 158 /xm lines, but also 
the CO, OH, H 2 0, CH + , and DCO + lines. Table 3 shows the 
setup for both line spectroscopy (henceforth linespec) and range 
spectroscopy (rangespec). The primary lines are highlighted 


in bold. The effective spectral resolution is 0.020 /x m at [0 1 ] 
63 /xm (88 km s -1 ) and 0.126 /tm for the [Cn] 158 /xm line 
(239 km s -1 ). Therefore the observed lines are unresolved in 
velocity even for the shortest wavelengths, which have the best 
spectral resolution, except in a few strong outflow sources, 
where the lines can be broad (Podio et al. 2012). 

For this study, spectra of 76 targets in 91 linespec observations 
were obtained, covering the [Oi] 63 /xm and DCO + 190 /xm 
lines (Table 4). The integration times for the majority of 
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Table 4 

GASPS Taurus Observations: PACS Spectroscopy 


Target 

R.A. (J2000) 
(h m s) 

Decl. (J2000) 
(deg m s) 

Obs. Mode 

Exposure 

(s) 

ObsID 

AA Tau 

4 34 55.420 

+24 28 53.16 

LineSpec 

1252 

1342190357 




LineSpec 

6628 

1342225758 




RangeSpec 

5141 

1342190356 




RangeSpec 

20555 

1342225759 

BP Tau 

4 19 15.830 

+29 06 26.90 

LineSpec 

1252 

1342192796 




LineSpec 

3316 

1342225728 

CIDA2 

4 15 05.157 

+28 08 46.21 

LineSpec 

1252 

1342216643 

Cl Tau 

4 33 52.000 

+22 50 30.20 

LineSpec 

1252 

1342192125 




RangeSpec 

5141 

1342192124 

CoKu Tau/4 

4 41 16.800 

+28 40 00.60 

LineSpec 

1252 

1342191360 




LineSpec 

6628 

1342225837 

CW Tau 

4 14 17.000 

+28 10 57.80 

LineSpec 

1252 

1342216221 




RangeSpec 

10279 

1342216222 

CX Tau 

4 14 47.860 

+26 48 11.01 

LineSpec 

3316 

1342225729 

CYTau 

4 17 33.730 

+28 20 46.90 

LineSpec 

1252 

1342192794 

DE Tau 

4 21 55.640 

+27 55 06.10 

LineSpec 

1252 

1342192797 




RangeSpec 

8316 

1342216648 

DF Tau 

4 27 03.080 

+25 42 23.30 

LineSpec 

1252 

1342190359 




RangeSpec 

5141 

1342190358 

DG Tau 

4 27 04.700 

+26 06 16.30 

LineSpec 

1252 

1342190382 




RangeSpec 

5141 

1342190383 

DG Tau B 

4 27 02.560 

+26 05 30.70 

LineSpec 

1252 

1342192798 




RangeSpec 

10279 

1342196652 

DHTau 

4 29 42.020 

+26 32 53.20 

LineSpec 

1252 

1342225734 

DI Tau a 

4 29 42.475 

+26 32 49.31 

LineSpec 

1252 

1342225734 

DK Tau 

4 30 44.240 

+26 01 24.80 

LineSpec 

1252 

1342192132 




LineSpec 

3316 

1342225732 




RangeSpec 

5141 

1342192133 

DL Tau 

4 33 39.060 

+25 20 38.23 

LineSpec 

1252 

1342190355 




LineSpec 

6628 

1342225800 




RangeSpec 

5141 

1342190354 

DM Tau 

4 33 48.720 

+18 10 10.00 

LineSpec 

1252 

1342192123 




LineSpec 

6628 

1342225825 




RangeSpec 

5141 

1342192122 

DN Tau 

4 35 27.370 

+24 14 58.90 

LineSpec 

1252 

1342192127 




LineSpec 

3316 

1342225757 




RangeSpec 

5141 

1342192126 

DO Tau 

4 38 28.580 

+26 10 49.44 

LineSpec 

1252 

1342190385 




RangeSpec 

5141 

1342190384 




RangeSpec 

12516 

1342225802 

DP Tau 

4 42 37.700 

+25 15 37.50 

LineSpec 

1252 

1342191362 




RangeSpec 

10279 

1342225827 

DQ Tau 

4 46 53.050 

+17 00 00.20 

LineSpec 

1252 

1342225806 




RangeSpec 

8316 

1342225807 

DS Tau 

4 47 48.110 

+29 25 14.45 

LineSpec 

3316 

1342225851 

FF Tau 

4 35 20.900 

+22 54 24.20 

LineSpec 

1252 

1342192802 

FM Tau 

4 14 13.580 

+28 12 49.20 

LineSpec 

1252 

1342216218 




RangeSpec 

10279 

1342216219 

FO Tau 

4 14 49.290 

+28 12 30.60 

LineSpec 

1252 

1342216645 




RangeSpec 

8316 

1342216644 

FQ Tau 

4 19 12.810 

+28 29 33.10 

LineSpec 

1252 

1342192795 




RangeSpec 

8316 

1342216650 

FS Tau 

4 22 02.180 

+26 57 30.50 

LineSpec 

1252 

1342192791 




RangeSpec 

10279 

1342194358 

FT Tau 

4 23 39.190 

+24 56 14.10 

LineSpec 

1252 

1342192790 

FW Tau 

4 29 29.710 

+26 16 53.20 

LineSpec 

1252 

1342225735 

FX Tau 

4 30 29.610 

+24 26 45.00 

LineSpec 

1252 

1342192800 

GG Tau 

4 32 30.350 

+17 31 40.60 

LineSpec 

1252 

1342192121 




RangeSpec 

5141 

1342192120 




RangeSpec 

12516 

1342225738 

GH Tau 

4 33 06.430 

+24 09 44.50 

LineSpec 

1252 

1342192801 




RangeSpec 

8316 

1342225762 
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Table 4 
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Target 

R.A. (J2000) 
(h m s) 

Decl. (J2000) 
(deg m s) 

Obs. Mode 

Exposure 

(s) 

ObsID 

GI Tau/GK Tau 

4 33 34.310 

+24 21 11.50 

LineSpec 

3316 

1342225760 

GM Aur 

4 55 10.990 

+30 21 59.25 

LineSpec 

1252 

1342191357 




RangeSpec 

5141 

1342191356 

GO Tau 

4 43 03.090 

+25 20 18.80 

LineSpec 

1252 

1342191361 




LineSpec 

3316 

1342225826 

Haro 6-5B b 

4 22 02.18 

+26 57 30.5 

LineSpec 

1252 

1342192791 

Haro 6-13 

4 32 15.410 

+24 28 59.70 

LineSpec 

1252 

1342192128 




RangeSpec 

5141 

1342192129 




RangeSpec 

12516 

1342225761 

Haro 6-37 

4 46 58.980 

+17 02 38.20 

LineSpec 

1252 

1342225805 

HBC 347 

3 29 38.370 

+24 30 38.00 

LineSpec 

1252 

1342192136 

HBC 356 

4 03 13.990 

+25 52 59.90 

LineSpec 

1252 

1342204134 




LineSpec 

1252 

1342214359 

HBC 358 

4 03 50.840 

+26 10 53.20 

LineSpec 

1252 

1342204347 




LineSpec 

1252 

1342214680 

HD 283572 

4 21 58.847 

+28 18 06.51 

LineSpec 

1660 

1342216646 

HK Tau 

4 31 50.570 

+24 24 18.07 

LineSpec 

3316 

1342225736 




RangeSpec 

8316 

1342225737 

HL Tau c 

4 31 38.437 

18 13 57.65 

LineSpec 

1252 

1342190351 

HN Tau 

4 33 39.350 

+17 51 52.37 

LineSpec 

3316 

1342225796 




RangeSpec 

10279 

1342225797 

HO Tau 

4 35 20.200 

+22 32 14.60 

LineSpec 

1252 

1342192803 

HV Tau 

4 38 35.280 

+26 10 38.63 

LineSpec 

3316 

1342225801 

IP Tau 

4 24 57.080 

+27 11 56.50 

LineSpec 

3316 

1342225756 

IQ Tau 

4 29 51.560 

+26 06 44.90 

LineSpec 

1252 

1342192135 




LineSpec 

3316 

1342225733 




RangeSpec 

5141 

1342192134 

IRAS 04158+2805 

4 18 58.140 

+28 12 23.50 

LineSpec 

1252 

1342192793 




RangeSpec 

5141 

1342192792 

IRAS 04385+2550 

4 41 38.820 

+25 56 26.75 

LineSpec 

3316 

1342225828 




RangeSpec 

8316 

1342225829 

Jl- 4872 

4 25 17.678 

+26 17 50.41 

LineSpec 

1660 

1342216653 

LkCal 

4 13 14.142 

+28 19 10.84 

LineSpec 

1660 

1342214679 

LkCa 3 

4 14 47.973 

+27 52 34.65 

LineSpec 

1660 

1342216220 

LkCa 4 

4 16 28.109 

+28 07 35.81 

LineSpec 

1660 

1342216642 

LkCa 5 

4 17 38.940 

+28 33 00.51 

LineSpec 

1660 

1342216641 

LkCa 7 

4 19 41.272 

+27 49 48.49 

LineSpec 

1660 

1342216649 

LkCa 15 

4 39 17.800 

+22 21 03.48 

LineSpec 

1252 

1342190387 




LineSpec 

6628 

1342225798 




RangeSpec 

5141 

1342190386 

RWAur 

5 07 49.540 

+30 24 05.07 

LineSpec 

1252 

1342191359 




RangeSpec 

5141 

1342191358 

RYTau 

4 21 57.400 

+28 26 35.54 

LineSpec 

1252 

1342190361 




RangeSpec 

5141 

1342190360 

SU Aur 

4 55 59.380 

+30 34 01.56 

LineSpec 

3316 

1342217844 




RangeSpec 

10279 

1342197845 

T Tau 

4 21 59.430 

+19 32 06.37 

LineSpec 

1252 

1342190353 




RangeSpec 

5141 

1342190352 




LineSpec 

1252 

1342215699 

UX Tau 

4 30 03.990 

+18 13 49.40 

LineSpec 

1252 

1342204350 




LineSpec 

1252 

1342214357 

UY Aur 

4 51 47.380 

+30 47 13.50 

LineSpec 

1252 

1342193206 




LineSpec 

1252 

1342215699 




RangeSpec 

10279 

1342226001 

UZ Tau 

4 32 42.890 

+25 52 32.60 

LineSpec 

1252 

1342192131 




RangeSpec 

5141 

1342192130 

V710 Tau 

4 31 57.800 

+18 21 35.10 

LineSpec 

1252 

1342192804 

V773 Tau 

4 14 12.920 

+28 12 12.45 

LineSpec 

3316 

1342216217 

V807 Tau d 

4 33 06.641 

+24 09 54.99 

LineSpec 

1252 

1342192801 




RangeSpec 

8316 

1342225762 

V819 Tau 

4 19 26.260 

+28 26 14.30 

LineSpec 

1660 

1342216651 

V836 Tau 

5 03 06.600 

+25 23 19.70 

LineSpec 

3316 

1342227634 
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Target 

R.A. (J2000) 
(h m s) 

Decl. (J2000) 
(deg m s) 

Obs. Mode 

Exposure 

(s) 

ObsID 

V927 Tau 

4 31 23.820 

+24 10 52.93 

LineSpec 

3316 

1342225763 

VI 096 Tau 

4 13 27.230 

+28 16 24.80 

LineSpec 

1252 

1342214678 

VY Tau 

4 39 17.410 

+22 47 53.40 

LineSpec 

1252 

1342192989 

XZ Tau 

4 31 39.480 

+18 13 55.70 

LineSpec 

1252 

1342190351 




RangeSpec 

5141 

1342190350 

ZZ Tau 

4 30 51.380 

+24 42 22.30 

LineSpec 

1252 

1342192799 


Notes. 

a DI Tau observed in the same field of view as DH Tau. 
b Haro 6-5B observed in the same field of view as FS Tau. 
c HL Tau observed in the same field of view as XZ Tau. 
d V807 Tau observed in the same field of view as GH Tau. 

observations was 1252 s; a few targets were re-observed with 
longer integrations ranging from 3316 s to 6628 s. Of the 
76 targets observed in linespec mode, 38 were also observed 
in rangespec mode (42 observations), with integrations of 
5141-10279 s, covering 70-200 \jl m, including the [O i] 145 /xm 
and [Cn] 158 /xm lines. The spectroscopic data were reduced 
with HIPE (Ott 2010) v4.2.0, utilizing the pipeline scripts 
provided at the Herschel data reduction workshop in 2010 
January. 18 

For each observation, we extracted the mean of the two nod 
positions to obtain flux values for each of the 25 spaxels per 
observation. We binned the spectra at half the instrumental 
spectral resolution (Nyquist sampling) for a given wavelength 
within HIPE, and output standard ASCII text files in order 
to extract continuum and line flux measurements. We utilized 
the standard IDL line fitting routine MPFITPEAK (Markwardt 
2009) to derive the total line flux and the continuum flux at the 
observed wavelength of each line. We fit a line to the observed 
continuum, defined nominally as the region of ~2-10 times the 
instrumental FWHM from the rest wavelength for each line. 
The error on the continuum is defined as the standard deviation 
of the residual about the mean. The reported line fluxes are 
the integrated flux of a Gaussian line fit, with errors calculated 
from the errors on the values of line amplitude and width, as 
reported by MPFITPEAK. For spectra with no line detection we 
report 3cr upper detection limits, where the lcr flux is calculated 
by integrating a Gaussian with height equal to the rms of 
the continuum and width equal to the instrumental FWHM. 
Typical rms error values for continuum levels in individual 
spaxel spectra are ~0.05 Jy. Typical line flux errors (rms) are 
~0.5-l x 10 -17 W m -2 . 

3.2.2. Correcting PACS Spectroscopy Observations 
for Pointing Errors 

The PACS instrument consists of an array of 25 “spaxels,” 
arranged 5 x 5 in a roughly square formation; see Poglitsch et al. 
Each spaxel subtends 9M x 9"4 on the sky. Ideally, observations 
are carried out with the celestial source being centered on the 
central spaxel. For a “well-pointed” observation, the largest flux 
measured in the spaxel array will, by construction, be located in 
the central spaxel. In such a case, the total flux of a point source 
can be estimated from the observed flux in the central spaxel 
by applying an aperture correction (which is the ratio of the 
observed flux in the spaxel to the total source flux, as estimated 


from the instrument PSF), as provided by the HIPE reduction 
software. 

However, some of the observations for our GASPS OT key 
program, especially for sources in Taurus, were affected by 
large errors in the pointing. In some observations the sources 
were found to be located up to 8" from the central spaxel, 
far greater than reported pointing accuracy 19 of ~2". As none 
of our Taurus sources are extended in continuum emission (cf. 
Section 3. 1), it is clear upon inspection of the relative continuum 
fluxes in the spaxels of each observation if a source is not well 
centered (Figure 1). In the case of the GASPS’ Taurus targets, the 
pointing solution shows the targets systematically to the east of 
the central spaxel. For these observations, the simple aperture 
correction method will not yield accurate results, especially 
since the reduction software provides no way of determining 
precisely where the source is located relative to the center of 
the spaxel array. Furthermore, the aperture correction formula 
makes use of the signal (information) in a single spaxel and 
ignores any observed flux in any additional spaxels. In order 
to determine accurately the flux of our targets from these 
mispointed observations and overcome the shortcomings of 
the aperture correction formula, we developed a method of 
estimating the offset of the source from the center of the spaxel 
array, and using that information to calculate an “optimum” 
value of the source flux. 

To do this, we simulated the observation of a source using 
wavelength-appropriate theoretical PACS PSFs 20 and the known 
spaxel sizes and locations relative to the central spaxel. The 
theoretical PSFs are very close to the observed PSFs (see 
Poglitsch et al. 2010), deviating from the observed PSFs on 
the order of only a few percent. Relative spaxel positions on the 
sky are obtained from the R.A. and decl. coordinates provided 
for each spaxel in HIPE. We take into account the position angle 
(PA.) of the telescope to transform the spaxel offsets from R.A. 
and decl. to arcseconds in (y, z) spacecraft coordinates. 

The simulations consisted of computing, for each spaxel, the 
fractional flux as well as the flux relative to the “observed” peak 
flux as a function of offset in spacecraft coordinates (y and z) 
between the peak of the PSF and the center of the spaxel array. 
We stepped the spaxel array across the wavelength appropriate 
theoretical PSFs, generating a library of spaxel fractional fluxes 
and relative fluxes for offsets of 3=20", in increments of 075, 
in both directions. For any given observation, the continuum 
fluxes in each spaxel were divided by the peak continuum value 


18 


https://nhscsci.ipac.caltech.edu/sc/index.php/~Workshops/ 
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1 9 http : //herschel . esac . esa. int/ twiki/bin/ view/Public/SummaryPointing 

20 http://pacs.ster.kuleuven.ac.be/pubtool.psf 
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Figure 1. “Footprint” for DG Tau. Spaxels are plotted as 1" on a side for clarity. Because DG Tau is not extended in continuum at 63 /zm, it is easily identified as falling 
off the central spaxel by examining the relative continuum values in each spaxel. Plotted as the green star is the best fit position of the target via y 1 minimization of 
spaxel continuum values and the theoretical PSF at 63 /zm. A error circle representing the error in the positional fit is plotted. The blue star is the SIMBAD position 
for DG Tau. 

(A color version of this figure is available in the online journal.) 


and then compared to the simulated relative fluxes at each offset 
position in the library using a/ 2 statistic. The positional offset 
(relative to the central spaxel) of the source is given by the offset 
location that yielded the minimum x 2 - As in the line fits, the 
region used to define the continuum in each spaxel was a range 
~2-10 times the instrumental FWHM from the rest wavelength 
for each line. 

Once the offset position was known, the fractional fluxes 
corresponding to that offset were used to estimate the optimal 
value of the total source flux (see, e.g., Horne 1986). For the 
best fit offset position, we calculated an optimal flux value via 


an error weighted average using only spaxels with a detected 
flux level above a signal-to-noise ratio (S/N) greater than 5. The 
optimal estimate is given by 


with an error of 


^opt. — 


E^O Wi-fi 

E "=0 W i 


^opt. — 




- 1/2 


(i) 


( 2 ) 
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The Wi are weights given by 


where pi is the fractional flux in spaxel i as given by the 
previous determination of the offset, and is the measured 
statistical error (the standard deviation on the mean continuum, 
as described above) on the continuum flux. For bright sources, 
several spaxels contribute to the final flux estimate, with each 
contribution weighted by the fractional flux in that spaxel for 
the estimated offset position. For weak sources, or those that 
are well centered in a single spaxel (and therefore, for which 
flux is detected in only a single spaxel), the method should give 
the same result as the simple aperture correction formula. A 
comparison of our “optimum” flux estimates and fluxes using 
the single spaxel and the aperture correction value provided 
by HIPE, for well-centered sources shows the values to be 
equivalent (i.e., for well-centered sources, or faint sources where 
flux is seen in only a single spaxel, the method of “optimal” 
flux estimates described herein is consistent with standard flux 
estimates utilizing only the HIPE reduction method). 

For targets in which two sources were present in the PACS 
field of view (ObsID 1342192193 containing FS Tau A/B, 
ObsID 1342190351 containing XZ and HL Tau, ObsID 
1342192801 containing V807 Tau and GH Tau), the method 
described above was modified somewhat to account for addi- 
tional flux from each source in each spaxel. The spaxels with 
the peak fluxes, corresponding to the approximate locations of 
the sources, were identified. For each individual source in the 
field of view, the method described above was carried out after 
setting the weights for all spaxels between the two peaks and 
those associated with the other source to zero. As above, the to- 
tal flux of the sources was then determined from a least- squares 
fit of the predicted fluxes in the remaining spaxels to the ob- 
served fluxes. This yielded accurate relative locations, and thus 
accurate continuum flux values, of the two sources in each field 
of view. 

Testing of the code was performed on sources falling off 
the central spaxel (as determined from “by-eye” inspection of 
relative continuum values) as well as on well-centered sources. 
The positional uncertainties in the y,z offsets, as defined by 
a lcr confidence contour based on the x 2 minimization, are 
less than 075 in radial offset, which is less than our step 
increments used in generating the fractional fluxes falling in 
each spaxel. As an additional check of the offset code, we 
find very good agreement in derived offset positions between 
different wavelengths (linespec and rangespec observations), 
which are separate observations, but usually done back to back 
and which show the same pointing anomalies. From a sample of 
19 targets obtained in both spectroscopy modes (line scan and 
range scan), we find agreement in derived radial offset positions 
computed using the 7 different wavelengths ranges, to be better 
than l" in radial offset. 

As described above, all spectroscopy continuum values re- 
ported herein are the result of the error weighted average of 
the aperture corrected flux in spaxels which have a S/N of 5 
or greater. Only one T Tauri star, T Tau itself, is extended at 
PACS wavelengths; millimeter interferometry of many of the 
Taurus sources (Andrews & Williams 2007; Isella et al. 2009; 
Guilloteau et al. 2011; Harris et al. 2012) show the disk di- 
ameters to be <2" on the sky. Since the observed sizes of the 
objects in our sample will be smaller at the PACS wavelengths 
of 60-200 /zm than at millimeter wavelengths, any “extended” 
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continuum emission in multiple spaxels above and beyond what 
would be expected from the PSF (i.e., adjacent spaxels with 
equal continuum measurements) is a result of mis-pointing of 
the telescope. 

However, the line emission can often be somewhat extended, 
especially for young T Tauri stars, as many of the sources in our 
sample are well-studied outflow sources. In disks surrounding 
stars with known outflows, summing the line flux contribution 
from all spaxels will overestimate the line emission arising from 
the disk, since the spectral and spatial resolution of PACS at 
these wavelengths is insufficient to resolve the outflow compo- 
nent of the line emission from that of the disk. Therefore, the 
observed [O i] 63 /zm emission in outflow sources only provides 
an upper limit to the emission from the disk. To minimize the 
contribution from the outflow, we only quote the line intensity 
from the central spaxel (well-centered observations) or the ex- 
pected line intensity at the continuum peak interpolated from 
the line intensities in nearby spaxels. This is only significant for 
spatially extended sources like T Tau and DG Tau B, where our 
line intensities are about a factor of two smaller than the total 
(i.e., spatially integrated) line intensities reported by Podio et al. 
( 2012 ). 

3.2.3. Comparison of Spectroscopy and Photometry 

We compare the continuum flux values obtained from pho- 
tometry (Table 1) with those from spectroscopy (Table 5) in 
Figure 2. The photometric calibration for PACS has been shown 
to have much higher absolute accuracy than the spectroscopic 
calibration (Vandenbussche et al. 2011; Muller et al. 2011). 
We find that the photometry continuum values are consistently 
higher by ~30-40% (rms ~ 20%) than the spectroscopic values 
at both 70 and 100 /zm. This calibration discrepancy appears 
to be wavelength dependent; comparison of the 160 /zm pho- 
tometric continuum and the 158 /zm spectroscopic continuum 
values show agreement to within 2% (rms ~ 20%), on average. 
We consider the absolute calibration of the spectroscopic data 
reduced with HIPE v4.2.0 to be accurate to within ~40%. 

In order to check how much the calibration has improved in 
the most recent HIPE release, vlO.O, we re-reduced the spectra 
at 63 /zm for a randomly selected sample of 15 stars ranging 
from very faint, i.e., continuum emission ^0.5 Jy at 63 /zm, 
intermediate (a few Jy), to bright (10 Jy) with about five stars 
in each category. The HIPE vlO.O calibration has certainly 
improved and is much closer to the results obtained from 
continuum photometry. From the comparison sample we find 
that the continuum flux densities are 18% brighter in HIPE vlO.O 
than compared to HIPE v4.2, while the line intensities are only 
11% brighter than in HIPE v4.2. Within errors, however, there 
is no clear difference between continuum and line calibration. 
For very faint sources HIPE vlO.O is no better than HIPE v4.2, 
for both releases the spectroscopy calibration can be off by 
more than 50% compared to each other and to very accurate 
photometry. 21 Since we have overlapping targets with the DIGIT 
project (Sturm et al. 2010), we also compared the targets in 
common; the well-pointed DIGIT observations agree extremely 
well with our photometry. In general we find agreement to better 
than 15% when comparing spectra reduced in HIPE v4.2 to that 


21 All spectroscopy and photometry will be re-reduced by the GASPS team as 
part of our data delivery once HIPE is stable. This is expected to change the 
overall calibration for the spectroscopy bands B2A and B2B (50-70 /xm) by 
~20%; for all other spectroscopy bands and for photometry the changes are 
expected to be small. 


11 


The Astrophysical Journal, 776:21 (25pp), 2013 October 10 


Howard et al. 


Table 5 

Continuum Flux in Spectroscopic Observations 


Name 

63.18/zm 

(Jy) 

72.84 /zm 

(Jy) 

78.74 /zm 

(Jy) 

90.16 /zm 
(Jy) 

145.53 /zm 
(Jy) 

157.74 /zm 

(Jy) 

179.53 /zm 
(Jy) 

189.57 /zm 
(Jy) 

AA Tau 

0.85 ± 0.02 

0.85 ± 0.02 

0.82 ± 0.02 

0.89 ± 0.01 

0.99 ± 0.01 

1.18 ±0.01 

1.12 ±0.02 

0.92 ± 0.05 

BP Tau 

0.45 ± 0.04 







<0.17 

Cl Tau 

1.38 ±0.05 

1.62 ±0.03 

1.57 ±0.05 

1.90 ±0.04 

2.04 ± 0.03 

2.37 ± 0.03 

2.52 ± 0.04 

2.06 ± 0.07 

CIDA2 

<0.09 







<0.24 

CoKu Tau-4 

0.90 ± 0.02 







0.61 ± 0.04 

CW Tau 

1.31 ±0.04 

1.48 ±0.03 

1.58 ±0.02 

1.59 ±0.03 

1.67 ±0.02 

1.93 ±0.02 

1.83 ±0.02 

1.81 ±0.10 

CX Tau 

0.30 ± 0.03 







<0.22 

CYTau 

<0.15 







<0.29 

DE Tau 

1.35 ±0.05 

0.70 ± 0.03 

0.67 ± 0.02 


0.65 ± 0.01 

0.62 ± 0.01 


<0.27 

DF Tau 

0.43 ± 0.05 

0.32 ± 0.03 

0.35 ± 0.04 

0.59 ± 0.04 

0.20 ± 0.02 

0.12 ±0.03 

0.20 ± 0.04 

<0.29 

DG Tau 

13.94 ±0.09 

11.94 ±0.04 

13.72 ± 0.07 

14.87 ± 0.04 

14.09 ± 0.03 

15.27 ± 0.03 

14.26 ± 0.05 

11.54 ±0.06 

DG Tau B 

9.85 ± 0.05 

10.37 ± 0.03 

10.96 ± 0.04 

11.79 ±0.03 

14.63 ± 0.02 

15.34 ± 0.02 

14.46 ± 0.03 

13.88 ± 0.09 

DHTau 

0.25 ± 0.05 







0.39 ± 0.07 

DITau 

<0.15 







<0.22 

DK Tau 

0.86 ± 0.03 

0.73 ± 0.04 

0.73 ± 0.04 

0.90 ± 0.05 

0.76 ± 0.03 

0.85 ± 0.03 

1.04 ±0.04 

<0.31 

DL Tau 

0.96 ± 0.02 

0.71 ± 0.03 

1.07 ±0.05 

1.08 ±0.05 

1.49 ±0.02 

1.77 ±0.03 

2.15 ±0.04 

1.58 ±0.03 

DM Tau 

0.58 ± 0.02 

0.68 ± 0.03 

0.86 ± 0.04 

0.81 ± 0.03 

0.81 ± 0.03 

0.93 ± 0.02 

1.04 ±0.04 

0.65 ± 0.03 

DN Tau 

0.69 ± 0.02 

0.55 ± 0.05 

0.70 ± 0.05 

0.58 ± 0.04 

0.63 ± 0.02 

0.77 ± 0.03 

0.78 ± 0.04 

0.70 ± 0.06 

DO Tau 

2.82 ± 0.05 

4.01 ± 0.02 

4.16 ±0.02 

3.15 ±0.04 

3.94 ± 0.01 

4.27 ± 0.01 

3.45 ± 0.04 

2.87 ± 0.07 

DP Tau 

0.79 ± 0.06 

0.36 ± 0.02 

0.32 ± 0.03 

0.37 ± 0.02 

0.18 ±0.02 

0.25 ± 0.02 

0.32 ± 0.03 

<0.31 

DQ Tau 

1.14 ±0.04 

1.13 ±0.02 

1.18 ±0.02 


0.98 ± 0.01 

1.09 ±0.02 


0.64 ± 0.09 

DS Tau 

0.17 ±0.03 







<0.17 

FF Tau 

<0.12 







<0.29 

FM Tau 

0.56 ± 0.06 

0.39 ± 0.02 

0.39 ± 0.02 

0.59 ± 0.03 

0.23 ± 0.02 

0.30 ± 0.02 

0.26 ± 0.02 

<0.29 

FO Tau 

0.42 ± 0.05 

0.35 ± 0.02 

0.39 ± 0.02 


0.23 ± 0.01 

0.26 ± 0.01 


<0.23 

FQ Tau 

<0.12 

0.11 ±0.02 

0.15 ±0.02 


0.06 ± 0.01 

0.09 ± 0.01 


<0.32 

FS Tau A 

2.34 ± 0.05 

1.73 ±0.03 

1.61 ±0.07 

1.88 ±0.03 

1.63 ±0.02 

1.78 ±0.02 

1.55 ±0.03 

1.47 ±0.09 

FT Tau 

0.52 ± 0.04 







1.06 ±0.11 

FW Tau 

<0.12 







<0.29 

FX Tau 

0.17 ±0.05 







<0.29 

GG Tau Aab 

2.93 ± 0.04 

3.97 ± 0.01 

4.51 ±0.03 

4.52 ± 0.04 

7.62 ± 0.01 

8.60 ± 0.01 

8.22 ± 0.03 

7.13 ±0.08 

GH Tau 

0.87 ± 0.05 

0.34 ± 0.02 

0.30 ± 0.02 


0.18 ±0.02 

0.20 ± 0.01 


0.31 ±0.07 

GI/GK Tau a 

0.91 ± 0.07 







0.82 ±0.12 

GM Aur 

1.99 ±0.05 

2.46 ± 0.05 

2.46 ± 0.04 

3.00 ± 0.04 

3.64 ± 0.02 

4.25 ± 0.03 

5.00 ± 0.04 

2.78 ± 0.07 

GO Tau 

0.24 ± 0.02 







0.50 ± 0.04 

Haro6-5B 

2.11 ±0.05 

1.87 ±0.02 

1.68 ±0.03 

1.79 ±0.03 

2.26 ± 0.02 

2.61 ± 0.02 

2.66 ± 0.03 

3.73 ±0.10 

Haro 6-13 

5.04 ± 0.05 

5.53 ± 0.05 

6.59 ± 0.02 

5.84 ± 0.04 

6.37 ± 0.01 

6.86 ± 0.02 

6.55 ± 0.05 

5.41 ± 0.08 

Haro 6-37 

0.85 ± 0.04 







1.42 ±0.07 

HBC 347 

<0.15 







<0.32 

HBC 356 

<0.09 







<0.24 

HBC 358 

<0.18 







<0.29 

HD 283572 

<0.09 







<0.22 

HK Tau 

2.03 ± 0.03 

2.13 ± 0.02 

2.29 ± 0.02 


2.13 ±0.01 

2.43 ± 0.01 


1.73 ±0.07 

HL Tau 

66.19 ±0.14 

65.15 ± 0.07 

65.98 ± 0.08 

67.79 ± 0.08 

65.91 ± 0.06 

69.37 ± 0.08 

53.90 ± 0.09 

45.45 ±0.11 

HN Tau 

0.77 ± 0.02 

0.71 ± 0.02 

0.79 ± 0.03 

0.74 ± 0.03 

0.60 ± 0.01 

0.64 ± 0.02 

0.60 ± 0.03 

0.55 ± 0.07 

HO Tau 

<0.12 







<0.34 

HV Tau C 

0.80 ± 0.04 







1.09 ±0.05 

IP Tau 

0.44 ± 0.03 







<0.16 

IQ Tau 

0.61 ± 0.03 

0.48 ± 0.06 

0.43 ± 0.03 

0.59 ± 0.04 

0.59 ± 0.02 

0.70 ± 0.03 

0.76 ± 0.04 

0.71 ± 0.06 

IRAS 04158+2805 

0.76 ± 0.05 

0.92 ± 0.03 

0.93 ± 0.04 

1.25 ±0.05 

2.45 ± 0.02 

3.18 ±0.03 

3.87 ± 0.03 

2.94 ± 0.08 

IRAS 04358+2550 

1.87 ±0.02 







1.44 ±0.04 

J 1-4827 

<0.09 







<0.27 

LkCal 

<0.09 







<0.26 

LkCa 3 

<0.12 







<0.27 

LkCa 4 

0.12 ±0.04 







<0.22 

LkCa 5 

<0.12 







<0.24 

LkCa 7 

<0.12 







<0.20 

LkCa 15 

0.99 ± 0.02 

0.98 ± 0.04 

1.08 ±0.04 

1.13 ±0.04 

1.55 ±0.03 

1.82 ±0.03 

2.18 ±0.06 

1.28 ±0.04 

RW Aur 

1.79 ±0.05 

1.85 ±0.04 

2.09 ± 0.04 

2.04 ± 0.04 

1.30 ±0.02 

1.54 ±0.03 

1.41 ±0.04 

1.12 ±0.13 

RYTau 

10.86 ± 0.07 

9.82 ± 0.03 

10.10 ±0.04 

10.00 ± 0.04 

7.98 ± 0.02 

8.64 ± 0.03 

8.50 ± 0.04 

5.73 ±0.11 

SU Aur 

5.95 ± 0.02 

5.44 ± 0.02 

4.72 ± 0.03 

4.11 ±0.03 

3.24 ± 0.02 

3.33 ± 0.02 

2.10 ±0.02 

1.49 ±0.05 

T Tau 

86.12 ±0.21 

77.48 ± 0.07 

75.93 ±0.10 

67.18 ±0.15 

39.53 ± 0.26 

39.64 ± 0.06 

35.77 ± 0.23 

25.28 ±0.13 

UX Tau 

3.28 ± 0.03 







2.31 ±0.07 

UY Aur 

4.85 ± 0.06 

4.87 ± 0.02 

4.89 ± 0.04 

4.52 ± 0.04 

3.11 ±0.01 

3.30 ± 0.02 

2.39 ± 0.02 

2.11 ±0.07 

UZ Tau 

1.23 ±0.06 

1.29 ±0.04 

1.47 ±0.05 

1.38 ±0.04 

1.78 ±0.03 

1.94 ±0.02 

1.69 ±0.04 

1.54 ±0.09 
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Table 5 

(Continued) 


Name 

63.18 /x m 

(Jy) 

72.84 /xm 

(Jy) 

78.74 /xm 

(Jy) 

90.16 /xm 
(Jy) 

145.53 /xm 
(Jy) 

157.74 /xm 
(Jy) 

179.53 /xm 

(Jy) 

189.57 /xm 
(Jy) 

V710 Tau 

0.47 ± 0.05 







1.10 ±0.18 

V773 Tau 

0.66 ± 0.03 







<0.17 

V807 Tau 

0.58 ± 0.05 

0.13 ±0.02 

0.12 ±0.02 


0.17 ±0.02 

0.18 ±0.02 


<0.29 

V819 Tau 

<0.12 







<0.22 

V836 Tau 

0.33 ± 0.02 







0.27 ± 0.04 

V927 Tau 

<0.09 







<0.17 

VI 096 Tau 

<0.15 







<0.27 

VY Tau 

<0.12 







<0.29 

XZ Tau 

4.20 ± 0.05 

3.98 ± 0.04 

4.04 ± 0.06 

3.89 ± 0.04 

2.90 ± 0.02 

2.78 ± 0.03 

2.87 ± 0.06 

3.71 ± 0.08 

ZZ Tau 

<0.12 







<0.29 


Note. a Gi Tau and GK Tau are unresolved in spectroscopy. 


of the DIGIT project, except for DG Tau, where the discrepancy 
is 50% (G. Herczeg 2011, private communication). 

The spectroscopic data in this paper are presented as a self- 
consistent dataset, completely reduced with HIPE v4.2 in order 
to avoid confusion regarding scaling values for targets reduced 
with different versions of HIPE. As the main results of this 
paper focus on comparison of continuum and line flux values 
at a given wavelength, the results should be unaffected by 
an overall scaling of line and continuum flux values, except 
possibly for line ratio values (discussed in Section 4.4). We note 
that in principle both the continuum and line intensities should 
be scaled by the appropriate value as determined by comparison 
to photometry for those wishing to use the continuum or 
line flux values contained herein. However, as there is not a 
100% overlap between sources observed in photometry and 
spectroscopy modes, and because of the wavelength dependence 
of the calibration, none of the spectroscopic continuum and line 
flux values (and corresponding plots) in this paper has been 
scaled to match those obtained by photometry. 

4. RESULTS OF SPECTROSCOPY 

Table 2 lists the [O i] 63 /z m, [O i] 145 /xm, and [C n] 158 /zm 
line flux values for our GASPS Taurus sample. The PACS 

spectroscopy continuum values are listed in Table 5. In some 

of our targets, other molecular lines like OH, CH + , and CO 

(for a complete listing, see Table 3) are also seen, primarily in 
known outflow sources (see, e.g., Podio et al. 2012). The 0 -H 2 O 
line at 63.323 /zm is also detected in eight of the sources in 

our sample (Riviere-Marichalar et al. 2012). This water line is 

only seen in outflow sources, although the emission is thought 
to originate in the disk rather than in the outflow. Although 

discussion of these molecular lines is outside the scope of this 
paper, targets which show evidence of these various molecular 
lines are noted in Table 2. Here we discuss only the strongest 
lines seen in the GASPS spectroscopic sample, i.e., the two [0 1 ] 
lines and the [C 11 ] line, which dominate the gas cooling (Tielens 

& Hollenbach 1985). 

4.1. [0 1 ] 63 fim Emission 

The [Oi] 63 /zm line was observed in 76 fields (Table 4) 

and detected in 43 stars at 3o or more. It was seen in all 

Class I objects, and in most Class II sources and transitional 

disk systems. We found that we always detected the [0 1 ] 63 /zm 
line at ^3cr if the 63 /z continuum flux was >0.5 Jy. The line 
is seen in known outflow sources as well as in “disk-only” and 


transitional disk systems. For outflow sources, the [0 1 ] 63 /zm 
line emission is essentially compact, whereas the jets seen in 
Ha or in forbidden optical lines often trace the outflow over 
one or several arcminutes from the star. The emission can be 
slightly extended in strong outflow sources (cf. Podio et al. 
2012). Line flux values in our sample range from 10“ 4 * * * * * * * * * * * * 17 W m -2 
(our sensitivity limit) to ~5 x 10“ 17 W m -2 for sources with no 
evidence of outflow and up to almost 10“ 14 W m -2 for known 
outflow sources. 

4.2. [0 1 ] 145 /zm Emission 

The [Oi] 145 /zm line, with an energy level of ~327 K, is 
more difficult to detect than the transition at 63 /zm. We detect 
the 145 /zm line at >3 a in 17 objects, nearly all of them outflow 
sources (Table 2); only one is a disk object. The non-detections 
of our disk sources show that the [0 1 ] 145 /zm line is more than 
10 times fainter than the [Oi] 63 /zm line in the disk. DE Tau 
is an exception, here the 145 /zm line is about as strong as the 
63 /zm line, see Section 5.7. This does not mean that the [Oi] 
145 /zm emission is absent in disks, but since our disk sources are 
relatively faint, the expected line emission is generally below our 
detection limit, see Section 4.4. Figure 3 shows three examples 
of the [0 1 ] 145 /zm line detections in our sample. 

4.3. [Cn] 158 /zm Emission 

For the [C 11 ] 158 /zm line the detection rate is even lower than 
for the [Oi] 145 /zm line; only 10 targets were detected, all of 
them outflow sources. Figure 4 shows three examples of [C 11 ] 
158 /zm spectra. Previous Infrared Space Observatory (ISO) 
observations (Creech-Eakman et al. 2002; Nisini et al. 2002; 
Liseau et al. 2006) show strong (~10 -16 to 10 -15 W m -2 ) [C 11 ] 
158 /zm emission in young (mostly Class 0 and Class I) stellar 
objects. The [Cn] is likely to originate in the outflow, rather 
than in the disk. In the sample of strong outflow sources studied 
by Podio et al. (2012), they found that the [Cn] emission was 
somewhat extended compared to the continuum emission from 
the disk for all but one source in their sample. In the prototypical 
outflow source HH 46 observed with PACS, van Kempen et al. 
(2010) found that the [C 11 ] emission was twice as strong in the 
blue- shifted outflow compared to the red- shifted outflow lobe 
or to the spaxel centered on the protostar, demonstrating that 
most, if not all of the [C 11 ] emission originates in the outflow 
rather than in the disk. 

For the outflow sources in our sample, we see [C 11 ] 158 /zm 
line fluxes of ~10 -17 W m -2 (DG Tau B and T Tau stand 
out with [Cn] 158 /zm line flux values of ~10 -16 W m -2 ), 
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Figure 2. Comparison of PACS photometric values obtained with HIPEv9 and 
PACS spectroscopic values obtained with HIPEv4.2. Top: 70 /xm photometry 
vs. 72.8 /xm spectroscopy. The photometric values at 70 /xm are systematically 
~42% (rms ~ 26%) higher than spectroscopy. Middle: at 100 /xm the photo- 
metric values are ~33% (rms ~ 20%) higher than spectroscopy. Bottom: at 
160 /xm and the photometric and spectroscopy values agree to within ~2% 
(rms ~ 20%), on average. 



144.5 145.0 145.5 146.0 146.5 

Wavelength (/xm) 


Figure 3. [O i] 145 line detections for three of the targets in the Taurus sample. 
The line at 144.784 /xm is CO / = 18 -> 17. 

(A color version of this figure is available in the online journal.) 



157.0 157.5 158.0 158.5 159.0 

Wavelength (/xm) 


Figure 4. [C n] 158 line detections for three of the targets in the Taurus sample. 
DQ Tau (bottom panel) is a non- detection. The line seen at <2cr is almost 
certainly spurious. 

(A color version of this figure is available in the online journal.) 

more than an order of magnitude lower than the previous ISO 
observations. Because the ISO beam was very large at 158 /xm, 
~80", it is likely that ISO picked up extended low surface 
brightness emission from the surrounding cloud or emission 
from unrelated sources. 


4.4. Line Ratios 

Figure 5 shows the line ratios for [Oi] 63 /xm, [Oi] 145 /xm, 
and [C n] 158 /xm. In all cases where the lines are detected at 3a 
or higher, the line ratios are 10-25 for [O i] 63/ [O i] 145. The line 
ratios for [Oi] 63/[Cn] 158 are similar. Due to the wavelength 
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Figure 5. [O i] 63 /zm vs. [O i] 145 /zm (left) and [C n] 158 /xm (right) line flux for the targets in the Taurus sample. Solid circles are sources for which there is a ^3<x 
detection for the lines. Sources for which the line detection is <3cr are plotted as 3cr upper limits. Red points are known outflow sources, black non-outflow, and green 
are transitional disks. The lower dashed line corresponds to a ratio of one, and the upper dashed line to a ratio of 10. Left: all of the outflow sources show ratios of 
~10 or higher, whereas the non outflow sources are typically undetected in [O i] 145 /zm, with the exception of DE Tau. 

(A color version of this figure is available in the online journal.) 


dependence on calibration with spectroscopy (discussed in 
Section 3.2.3), these ratios are likely to be even higher, as 
the [O i] 63 /an line flux is systematically underestimated in 
comparison to the [O i] 145 /xm and [C ii] 158 /im line fluxes. 

4.5. 63 /xm Continuum versus [Oi] 63 fim Emission 

Since the [O i] 63 /xm line can be strong in outflows, we 
divided our sample into two sets: (1) the “disk-only” stars, 
where the star has no or only weak outflow activity and (2) 
jet/outflow sources, where there is clear evidence for outflow 
activity (see Section 2). In Figure 6 we plot the [O i] 63 /xm 
line intensities versus the 63 /xm continuum flux density. We 
find a tight correlation between [O i] line emission and 63 /xm 
continuum flux with a correlation coefficient of 0.90 for the 
disk sources which have both line and continuum fluxes >3cr. 
A weighted fit gives 

log(Oi) = (0.737 ± 0.06) x log(S 63 ) - (0.22 ± 0.02), (4) 

where Oi is the [Oi] 63 /xm line flux in 10“ 16 W m -2 and 
S (53 is the continuum flux in Jy at 63 /xm. This fit is shown 
as a black line in Figure 6. Our typical 3 a line detection 
limits, ~10 -17 W m -2 , are not stringent enough to conclusively 
determine if the trend continues to values below ~0.3 Jy. All 
transitional disks lie below the correlation. In order to verify this 
correlation is valid beyond the 12 detected disk-only sources, we 
have used Kendall’s Rank Correlation, a correlation statistic that 
accounts for non-detections, but does not include uncertainties 


on detected sources (Isobe et al. 1986). We use an IDL adaptation 
of the algorithm presented in Isobe et al. (1986) to carry out 
the calculation. The population of 12 detected sources have a 
Z- value of 3.57, indicating there is a 3.6 x 10 -4 probability 
that a correlation is not present. When all 40 disk-only objects 
are included, the Z- value increases to 5.29, indicating that the 
probability of no correlation is essentially 0. 

Furthermore, we carry out a log-log fit to the fluxes and 
3 a upper limits of all disk-only sources. This fit is a modified 
X 2 minimization using the Levenburg-Marquardt algorithm as 
implemented in the IDL code MPFIT (Markwardt 2009). For 
detections, the contribution to x 2 is determined as normal (i.e., 
((O — M)/cr) 2 ). For non-detections, the x 2 value is set such that 
if the model value is below the 4cr limit, the point contributes a 
value of 1 to the deviates. For values above 4cr, the x 2 value is 
calculated as if the objected were detected at the 3a value with an 
uncertainty of lcr. Parameter uncertainties from this procedure 
are lower-limits, and do not reflect the two-sided nature of the 
uncertainties. 

This log-log fit results in a slightly shallower slope, 0.69 ± 
0.04, but is still well within the lcr errors. These tests of the 
correlation and the line-fit to the disk-only locus confirm that 
the line flux and the 63 /xm flux are strongly correlated. 

While the outflow sources show a similar trend, i.e., that 
sources with higher continuum flux tend to have higher line 
flux, they can exhibit line fluxes that are up to ~20 times higher 
than non-outflow sources for a given continuum value. For these 
outflow sources, it is clear that the [O i] 63 /xm line emission is 
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Figure 6. [Oi] 63 /zm line emission vs. 63 /xm continuum emission for the 
GASPS Taurus sample. Non-outflow sources are plotted in black, outflow 
sources are plotted in red, and transitional disks are plotted in green. Sources 
for which there is less than 3cr detection in line flux or continuum are plotted 
in gray as 3 o upper limits. None of the sources plotted in gray are known to 
drive outflows. There is a tight correlation (lower black line) between [O i] line 
emission and 63 /zm continuum flux for non-outflow sources (sources with less 
than 3 o detection in line or continuum flux (the known outflow sources and the 
transitional disks were omitted from the line fit), suggesting that the emission 
originates from the same part of the disk. The [O i] 63 /zm line emission in 
outflow sources is dominated by the outflow, and can be up to 20 times stronger 
(upper dashed line) than the emission from the disk. 

(A color version of this figure is available in the online journal.) 


completely dominated by the outflow. However, there are some 
outflow sources, which lie along the best-fit line to the disk-only 
sources (6 out of 24), suggesting that in some cases the shocks 
in the jets are not strong enough to excite the [O i] 63 /xm. 
Nevertheless, the strong correlation between [O i] emission and 
63 /xm continuum for disk sources suggests that one may can use 
low spectral resolution [O i] observations to determine whether 
a PMS star drives an outflow or use the derived correlation to 
provide a rough estimate of how much of the [O i] 63 /xm line 
emission originates in the disk. 

Several sources in our sample were detected in continuum 
emission at 63 /xm but showed little (below 3a) or no [Oi] 
63 /xm line emission in our first set of observations. Some of the 
sources for which there was a hint of a line were re-observed 
with roughly three times the exposure, yielding line detections 
on the order of 0.5-1 x 10“ 17 W m -2 as shown in Figure 6. 
All the sources, which were re-observed, were detected. For 
several re-observed sources, the line flux is lower than what 
the correlation would predict. Some of these sources, however, 
are known transitional systems (DM Tau, DN Tau, UX Tau 
A, GM Aur, and LkCa 15), which we already identified as 
having weaker [Oi] line intensities than CTTSs for the same 
continuum flux. This result suggests that the correlation between 
line flux and continuum is not only a good indicator of whether 
a source drives an outflow, but can also serve as a diagnostic 
for transitional disks. These transitional disks will be discussed 
in detail in a forthcoming paper (F. Menard et al. 2013, in 
preparation). 


If we include marginal line detections, i.e., >2a but <3 a, we 
find an additional seven sources which lie below the correlation 
for “disk-only” sources. Four of these marginal detections are 
transitional disks (FO Tau, CX Tau, IP Tau, and V836 Tau; 
Najita et al. 2007; Espaillat et al. 2011; Furlan et al. 2011), 
which confirm the trend we already saw, i.e., transitional disks 
are intrinsically fainter in [O i] 63 /xm than classical disks. The 
remaining three (Haro 6-37, FM Tau, and V7 10 Tau; see Table 2) 
may be entering the transitional phase of their evolution, or 
alternatively they appear faint due to increased systematic errors 
for faint sources. 

One other source stand out, GH Tau. It has a continuum flux 
density of 0.8 Jy at 63 /xm (Table 5), yet the [Oi] 63 /xm line 
emission is less than 3a. However, at 70 /xm the flux density 
is only 0.37 Jy (Table 1). Therefore the continuum emission at 
63 /xm is overestimated by more than a factor of two, and the 
[O i] 63 /xm 2a detection is at the level we would expect. 

5. DISCUSSION 
5.1. [O i] 63 jim Emission 

The [O i] 63 /xm line, with an upper energy level of 228 K, is 
one of the strongest FIR lines observed in the PACS wavelength 
range and a dominant cooling transition. Emission is thought to 
arise from the surface region of almost the entire disk (Gorti et al. 
2011; Kamp et al. 2011) with 50% of the emission coming from 
outside of 100 AU (Kamp et al. 2010). Meijerink et al. (2012) 
find the bulk of the atomic oxygen to be in LTE, and hence the 
[O i] 63 /xm and 145 /xm lines should be sensitive to the average 
temperature of the emission regions they are probing, making 
them potential probes of the energy of the gas. Several models 
show that the [O i] 63 /xm line is optically thick everywhere in 
the disk (Gorti & Hollenbach 2008; Kamp et al. 2011). 

Thermo-chemical models of disks (Gorti & Hollenbach 2008; 
Kamp et al. 2011) predict much higher line fluxes than what 
we observe for our disk sources, even after accounting for 
calibration uncertainties (our flux estimates are predicted to be 
~30% low, see Section 3.2.3), suggesting that the [O i] emission 
is still not fully understood. Outflow sources have much stronger 
line fluxes than their disk-only counterparts, with values that can 
be up to 20 times stronger for the same continuum flux compared 
to a non-outflow source. Therefore an outflow, if present, can 
dominate the line emission, especially for strong outflows. This 
is discussed further in Section 4.5. 

5.2. What Can We Learn from the [Oi 63/145 Line Ratios? 

If the emission of [Oi] 145 /xm and [Oi] 63 /xm came from 
the same region, our observed line rations 10-25, would suggest 
optically thin emission originating from gas with a temperature 
of a few 100 K (Tielens & Hollenbach 1985) located at the inner 
part and surface layers of the disk or from the shock regions in the 
outflow. However, we would also see similar ratios for optically 
thick lines, if the [Oi] 145 /xm line, which requires hotter gas, 
came from a more compact source (like only the inner part of 
the disk). 

Disk sources (DE Tau is an exception) and transitional disks 
do not show either [Oi] 145 /xm or [Cn] 158 /xm, because 
these disks are intrinsically faint and the emission is below 
our detection limit. 

Kamp et al. (2011) explored the [Oi] 63/145 line ratios 
using the results from the DENT grid (Disk Evolution with 
Neat Theory) that consists of 300,000 disk models with 1 1 free 
parameters. They find that the median of the [Oi] 63/145 line 
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Figure 7. Disk mass estimates from Andrews & Williams (2005) as a function of [O i] 63 /zm line flux and 63 /zm spectroscopic continuum flux (right panel) and 
the 63 /zm continuum (left panel). Detections of <3cr in either continuum or line flux are omitted. Red points are targets with known outflow activity, black points 
are sources with no known outflow, and green points represent known transitional disk sources. Known outflow sources have been omitted from line flux plot, as our 
analysis show that when a star drives an outflow, the outflow dominates the [O i] 63 /zm emission. There is no correlation between disk mass and [O i] 63 /zm line 
emission nor with the 63 /zm continuum flux. 

(A color version of this figure is available in the online journal.) 



ratios for the whole DENT grid is 25, while the ratio is ~16 
for the canonical gas-to-dust ratio of 100. Similar ratios are 
found for the disk models by Meijerink et al. (2008), who found 
that [O i] 145 /zm line emission is 20-40 times fainter than the 
[O i] 63 /z m line. This also agrees with the results by Meeus 
et al. (2012), who studied a sample of isolated HAEBE stars. 
They found that when the [O i] 145 /tm line was detected, it was 
typically 20-30 times fainter than the [O i] 63 /xm line. For such 
line ratios we would not expect to detect the [Oi] 145 /xm line 
in any of our “disk-only” sources. Line ratios of 10-50 are also 
predicted by shock models (Hollenbach & McKee 1989; Flower 
& Pineaut des Forets 2010). Thus the [Oi] 63/145 line ratios 
does not help us to discriminate between a jet or a disk origin. 

Liseau et al. (2006), who analyzed all ISO - LWS observations 
of young objects with outflows (mostly Class I and Class 0 
objects), found surprisingly low [O i] 63 /xm to [O i] 145 /xm 
ratios; more than half of their sample had line ratios less than 
10. The average line ratio of the ones that had a line ratio less than 
10 was 5.3 zb 2.6. Even though our sample consists mostly of 
Class II sources, we do include some Class I sources and several 
of our Class II sources drive powerful outflows. We would 
therefore have expected to see some low ratios, especially since 
all the sources detected in [O i] 145 /xm are outflow sources, but 
we do not. Essentially all of our detections have a ratio >10, see 
Figure 5. A recent paper from the guaranteed time Key program 
WISH (Karska et al. 2013) discuss PACS observations of [Oi] 
in 18 low-mass protostars, 5 of which were Class I and the rest 
were Class 0 sources. They also find much higher line ratios 
than Liseau et al. (2006). In their Class 0 sample Karska et al. 
found the [Oi] 63 /xm to [Oi] 145 /xm line ratios to vary from 
5.5 to >45, with a median of ~10.4, while their Class I samples 
varies from 10.4 to 26.5. Overall there does not seem to be much 
of a difference between the results obtained by Karska et al. and 
us. Their Class 0 sample has somewhat smaller ratios but the 
ratios for the Class I sources are about the same as we find in 


Taurus. The most likely explanation for the lower ratios found 
by Liseau et al. is probably due to calibration. The difference 
in beam size should have a rather minor effect, since the [O i] 
emission is rather compact, and if anything one would expect 
the 63 /x line to be spatially more extended than the 145 /xm line, 
which would push the ratio to higher values, not lower ones. 

5.3. Disk Mass versus [Oi] 63 fim Line Emission 

Kamp et al. (2010) suggested that the [Oi] 63 /xm and the 
[Cn] 158 /xm lines, especially when combined with CO, may 
provide a tool to measure the disk gas mass (see also, Kamp et al. 
2011). Our observations, however, show that even with Herschel 
we do not have the sensitivity to detect [C n] in disks around 
low-mass stars; all our [Cn] detections are outflow sources 
and the line emission, when detected, is almost certainly from 
shock or PDR emission in the outflow (Podio et al. 2012). On 
the other hand, [O i] 63 /xm, as shown in the previous section, 
clearly originates in the disk, because we detected the line in 
20 stars, which show little or no outflow activity. There are 
not enough observations of any given CO transition to enable 
a meaningful investigation of how the [Oi]/CO ratio correlates 
with disk mass. However, if the [O i] 63 /xm emission probes 
the disk gas mass, one might expect some sort of correlation 
between [O i] line flux and total disk mass. 

After searching through the literature, we decided to use 
the disk mass estimates for the large sample investigated by 
Andrews & Williams (2005). A homogeneous treatment is 
more important than accuracy for this investigation, because 
published mass estimates vary greatly depending on how they 
have been estimated, and what dust emissivities have been used 
(see, e.g., Ricci et al. 2010, 2012; Guilloteau et al. 2011). As we 
can see from the right panel in Figure 7, there is no correlation 
between [Oi] 63 /xm line flux and disk mass. We increased 
our sample by plotting the disk mass as a function of 63 /xm 
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continuum flux density, since the continuum emission originates 
in the disk and is not affected by the outflow. However, as we 
can see from the left panel in Figure 7, where we plotted outflow 
sources in red, they also reveal no correlation and the same is 
true for the whole sample. For the same 63 /x m continuum flux 
density or line intensity the disk mass can vary by more than a 
factor of 100. 

5.4. Source of 63 jim Continuum Emission 

Why is there a correlation between [O i] 63 /xm and the 63 /xm 
continuum? The simplest explanation for the tight correlation 
is that both line and continuum emission come from the same 
region of the disk. In order to investigate where the 63 /xm 
continuum emission originates in the disk, we first investigated 
the simple disk model used by Andrews & Williams (2007), 
which uses the broadband SEDs and sub-millimeter visibilities 
to constrain some of the basic parameters that describe the 
structure of the disk. Using their flat disk model (as described 
in Andrews & Williams 2007), we find that for the majority 
of their sources that are also in our Taurus sample, 90% of 
the 63 /xm continuum emission comes from the inner ~10%, 
or ~5-50 AU, depending on the size of the disk. It should be 
noted, however, that overall disk size is poorly constrained in 
these models. Models of T Tauri stars with MFOST (Pinte et al. 
2006), which include disk flaring, give similar results. For the 
Class II objects in our Taurus sample, it therefore appears that 
a significant fraction of the line emission also originates in the 
inner part of the disk. 

As a next step we explored the bulk properties of the 
disks using isothermal two-component graybody models. For 
these fits we use all published millimeter and sub-millimeter 
continuum fluxes combined with the far-infrared fluxes from 
this paper. To constrain the warm dust in the disk we also 
included MIPS 24 /xm and/or WISE 22 /xm data from the WISE 
All- Sky survey. The MIPS 24 /xm data are taken from (Luhman 
et al. 2010) and/or Rebull et al. (2010); the latter also includes 
MIPS 70 /xm data. The 22 /xm and 24 /xm mid-IR flux densities 
are still completely dominated by dust emission from the disk 
for the Taurus Class II sources and the contribution from 
the photospheric emission is negligible. The fitting program 
allows us to put in constraints for the dust temperature, and fits 
separately the size of both the cold and the warm dust as well 
as the dust emissivity index. In order to limit the number of 
fitted parameters we do not fit the dust emissivity index, f , of 
the warm component; instead we set it to 1. Whether we use 

= 1.0 or 1.5 does not really matter, the dust temperature is 
essentially the same to within a degree. The results of these fits 
for a sample of 21 sources (both disk and outflow sources, but not 
transitional disks), for which we had sufficient data, indicate that 
the temperature of the warm dust ranges from ~95 K to 190 K, 
with a median of 130 K. This is about the temperature we 
would expect for the [Oi] emitting gas (U. Gorti 2012, private 
communication). There is no difference in the temperature of the 
warm dust between outflow sources and non-outflow sources. 
Figure 8 presents several examples of these graybody fits. These 
results suggest that the [O i] emission originates in regions where 
the gas and dust are approximately thermalized. This might 
be the inner part of the disk or the warm, lower surface layers 
of the disk. 

As we saw earlier, we also see a correlation between [Oi] 
63 /xm flux and continuum emission in the transitional disks, 
but the [Oi] emission is consistently weaker by a factor of 
two or more for the same continuum flux. Two component 





Figure 8. Representative two component graybody fits. Millimeter and sub- 
millimeter data come from literature, far-infrared data from this paper, MIPS 24 
and/or 70 /xm fluxes are also from literature, while WISE 22 /xm fluxes are from 
the WISE All-SKy Catalog. Since we do not include photometry shortward of 
12 /xm, we can ignore the hot inner disk and the stellar photosphere. 

graybody fits of these transitional disks yield substantially cooler 
temperatures for the warm dust, i.e., the dust dominating the 
63 /xm continuum emission. For this sample the warm dust 
covers a relatively narrow temperature range, about 80-90 K, 
with a median of 85 K. This is consistent with the lower observed 
[O i] 63 /xm line emission in these transitional disks. If the dust 
and gas is coupled in [Oi] emitting regions of the disk, the 
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Figure 9. Cumulative distribution of separations for binaries systems in our 
sample. The black histogram represents all multiple systems in our survey. 
Dashed and dotted histograms represent subsamples of object which are 
respectively detected and non-detected in 70 /xm (red) and 850 /xm (green) 
continuum and [Oi] 63 /xm (blue). Tight binaries (below ~50 AU) are more 
likely to be undetected in all three types of observations than wider systems. 

(A color version of this figure is available in the online journal.) 


lower temperature compared to “normal” Class II disks readily 
explains the lower line intensity. We note, however, that this 
interpretation may be oversimplified. More detailed modeling 
of transition disks will be presented in a forthcoming paper 
(F. Menard et al. 2013, in preparation) showing that lowering 
the dust temperature alone makes the line too low (i.e., does 
not change the ratio properly) and that direct illumination of the 
disk rim probably needs to be taken carefully into account as 
well. 


5.5. Multiplicity 

The presence of stellar companions has potentially dramatic 
consequences on the properties of circumstellar disks through 
the tidal forces exerted on them. For instance, it has been 
established that T Tauri stars in close binaries (<50 AU) are 
much less likely to host disks (Cieza et al. 2009; Kraus et al. 
2012). Furthermore, there is evidence that disks in these systems 
are much less massive than those in wider binaries or around 
single stars (Harris et al. 2012, but see Duchene 2010). It is 
thus natural to explore the connection between multiplicity 
and our continuum and line observations. In particular, our 
survey allows one to address this question for the gaseous 
component of protoplanetary disks. To this end, we have 
compiled the multiplicity properties of each system in our 
sample (Kenyon et al. 2008; Kraus et al. 2011; R. J. White 
et al. 2013, in preparation). For high-order multiple systems, 
we first evaluated which component hosts the disk dominating 
the far-infrared emission, typically from high-resolution sub- 
millimeter mapping (e.g., Harris et al. 2012), and assigned the 
separation of the system to the closest companion around that 
component but excluding spectroscopic companions. 

The detection rate of multiple systems is lower than that of 
single stars for both continuum and line emission. The con- 
tinuum detection rate is 32/56 for binaries and 36/55 for sin- 
gles with a significance level of 90%. the line detection rate 
is 20/41 for binaries and 23/34 for singles with a significance 
level of 95%. Not surprisingly, we find that systems tighter than 
50-100 AU are less likely to be detected at far-infrared con- 
tinuum or [Oi] 63 /im emission than wider systems (Figure 9), 



Figure 10. Far-infrared and sub-millimeter continuum and line fluxes for all 
multiple systems in our sample as a function of the system’s (projected) 
separation. For high-order multiple systems, we assigned to the system the 
separation to the closest companion to the disk-bearing component of the system 
to the exclusion of spectroscopic companions. We plot here the continuum 
fluxes at 70 /xm (PACS, red diamonds), 160 /xm (PACS, purple squares), and 
850 /xm (Andrews & Williams 2005, green crosses), and the [Oi] 63 /xm line 
flux (PACS, blue triangles). Uncertainties displayed on PACS datapoint do not 
include calibration uncertainties, whereas a typical uncertainty of 10%-20% 
applies to all sub-millimeter points. Notice how binaries tighter than 50-100 AU 
have systematically reduced sub-millimeter fluxes whereas their far-infrared 
continuum and line fluxes span the same range as those of wider systems. 
The “error bars” on the right-hand side of the plot show the median and the 
34 percentile on each side (e.g., lex range) for single stars in our sample. 

(A color version of this figure is available in the online journal.) 


similar to conclusions reached from near- and mid-infrared as 
well as sub-millimeter wavelengths. Taken at face value, this 
suggests that the gas and dust components of disks are similarly 
affected by the presence of a companion. We note, however, that 
the detailed interpretation of these conclusions is complicated by 
the biased nature of the GASPS sample, which over-represents 
disk-bearing systems, especially for spectroscopic observations. 
Nonetheless, an intriguing finding is that neither the 70 /xm nor 
the [O i] 63 /xm line flux show a correlation with the binary 
separation in systems for which they are detected (Figure 10). 
This is opposite the results of sub-millimeter surveys (e.g., 
Harris et al. 2012) but in line with near- and mid-infrared studies 
(Cieza et al. 2009), further suggesting that both the continuum 
and line emission arise from the inner most regions of the disk, 
which are not truncated by stellar companions located 10 AU or 
more from the disk-bearing star. 

5.(5. Class III Objects (Weak-line T Tauri Stars) 

We do not detect line emission in any of the 15 WTTSs, 
which we observed in [Oi] 63 /xm (Table 2). In continuum we 
observed 44 Class III objects and detected only three of them: 
FW Tau, V397 Aur, and V819 Tau. All three have 70 /xm flux 
densities far below our [O i] 63 /xm detection threshold, if they 
follow the correlation found for “disk-only” sources, and would 
therefore require unrealistically long integration times to reach 
the expected signal level we predict from CTTS. V819 Tau 
(spectral type K7) was classified as Class II by Luhman et al. 
(2010) and Class III by Rebull et al. (2010). It was detected at 
70 /xm by Cieza et al. (2013), who did detailed SED modeling 
showing that the star has no excess in the IRAC bands but a 
clear excess at 24 /xm and 70 /xm. Their modeling suggests that 
V8 1 9 Tau has a warm debris disk. We detected V8 1 9 Tau at both 
70 /xm and 100 /xm, because we have longer exposure times than 
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Right Ascension (2000.0) 

Figure 11. False color images of CoKu Tau/4 overlaid with logarithmic contours. CoKu Tau/4 is surrounded by an extensive far-infrared nebulosity at 160 /im, which 
is much fainter at 100 /xm. 

(A color version of this figure is available in the online journal.) 


Cieza et al. We see no trace of emission at 160 /tm. Both FW 
Tau and V397 Tau were detected in all bands observed. Whether 
these disks are at the end of their primordial phase or whether 
the stars have already reached the debris disk stage is impossible 
to judge without detailed modeling. 

5. 7. Comments on Individual Objects 

Here we highlight a few stars, which previously were not 
known to power jets or which appear to have unusual character- 
istics. 

AA Tau, spectral type K7-M0 V, is a well-studied CTTS, 
which has long been believed to power an optical jet (Hirth 
et al. 1997; Bouvier et al. 1999, 2003, 2007). This was recently 
confirmed by Cox et al. (2013) who found a poorly collimated, 
faint jet extending to 21" from the star. The [Oi] 76300 line in 
AA Tau is narrow and slightly blue- shifted (Hirth et al. 1997; 
Hartigan et al. 1995). The Hei 710830 line is likewise faint, 
narrow and slightly blue-shifted (Edwards et al. 2006), whereas 
it is expected to be deep, broad and blue- shifted for an outflow 
source (cf. discussion on DL Tau below). 

Baldovin-Saveedra et al. (201 1) searched for mid-IR emission 
lines in a sample of 64 PMS stars in Taurus and detected [Ne ii] 
12.81 /zm emission in 18 objects including AA Tau. They found 
in general that the luminosity of the [Neu] line was stronger 
for sources driving jets than that for those without known jets, 
although not for AA Tau. Najita et al. (2009), who analyzed 
high-resolution [Ne n] emission line profiles from AA Tau found 
that they could be explained by originating in the disk, although 
they could not rule out contribution from a jet. Our [O i] 63 /zm 
observations, however, show that [Oi] line is consistent with 
a disk origin (Figure 6), i.e., with little or no excess emission 
from the jet. However, since it has been shown to drive a faint 
jet, we classify it as a jet/outflow source. The inner disk of 
AA Tau is rich in organic molecules (Carr & Najita 2008). The 
0 -H 2 O line at 63.32 /zm has also been detected for this source 
(Riviere-Marichalar et al. 2012). 

CoKu Tau/4 is a Ml. 5 T Tauri star with weak Ha emission 
(Cohen & Kuhi 1979; Kenyon et al. 1998) and it therefore has 
a low accretion rate. It was identified by Forrest et al. (2004) 
as a transitional disk based on Spitzer / IRS observations, which 
revealed that it has no infrared excess shortward of 8 /zm, but 
a large excess at 20-30 /zm indicating that the region within 
10 AU have been largely cleared of dust. Ireland & Kraus (2008) 
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Figure 12. WFPC 2 image (0.8 /xm) of the reflection nebulosity illuminated 
by CoKu Tau/4. The PACS 160 /xm image (contours) closely follows the 
nebulosity. 

(A color version of this figure is available in the online journal.) 

carried out aperture-masking interferometry and adaptive optics 
imaging of the star and show that CoKu Tau/4 is a near-equal 
binary star of projected separation ~8 AU (53 mas). The binary 
has therefore cleared the inner disk. Both stars are of spectral 
type M1-M2 and are ~3-4 Myr old (Ireland & Kraus 2008). It 
is therefore not a transitional disk. 

CoKu Tau/4 is a strong far-infrared source (Table 1) and 
one of the few stars in our sample which is associated with a 
far-infrared nebulosity (Figure 11). At 70 /zm the star appears 
point-like, but at 100 /zm nebulosity can be seen to the south 
and southeast of the star. At 160 /zm the emission extends over 
more than an arcminute (Figure 11). The far-infrared nebulosity 
follows closely the reflection nebulosity seen with WFPC 2 
camera on the Hubble Space Telescope 22 (Figure 12). The [Oi] 
flux is consistent with a disk origin (Figure 6). 

DE Tau is a young CTTS of spectral type M2 (Kenyon & 
Hartmann 1995). It is not known to drive an outflow. It stands out 


22 The WFPC 2 image for CoKu Tau 4 was obtained as part of HST program 
GO 9160; PI: D. L. Padgett. 
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in our sample as the only non-outflow star with a 3a detection 
of the [Oi] 145 /x m line and only a marginal (~2 a) detection 
of the [O i] 63 /xm line. For DE Tau the upper limit of the 
line ratio of [O i] 63 to [O i] 145 /xm is on the order of unity, 
suggesting that both transitions are optically thick and originate 
from the same part of the disk. An isothermal graybody fit to the 
millimeter and far-infrared SED requires the dust to be optically 
thick at millimeter wavelengths, i.e., a dust emissivity index of 
0, suggesting that DE Tau could be surrounded by a compact 
optically thick disk. 

DG Tau and DG Tau B both power optical jets (Mundt & 
Fried 1983; Mundt et al. 1987; Kepner et al. 1993; Eisloffel & 
Mundt 1998). DG Tau is a young, heavily accreting star with 
a spectral type of K7-M0 (Gullbring et al. 1998), obscured by 
a visual extinction of 3.2-5 .4 mag (Fischer et al. 2011). It also 
powers a low velocity, wide-angled bipolar molecular outflow 
(G. Sandell, 2013, private communication). The disk has been 
imaged with high spatial resolution at millimeter wavelengths 
and is seen relatively face on with an inclination of ~30° (Isella 
et al. 2010; Guilloteau et al. 2011). In this case we find the [O i] 
63 /xm emission to be completely dominated by the outflow 
but only marginally extended. Podio et al. (2013) spectrally 
resolved the [Cn] 158 /xm line with HIFI, which shows that 
the [C n] emission originates in the blue-shifted jet, and not in 
the disk. 

DG Tau B is a deeply embedded Class I source. The spectral 
type is unknown, but modeling of the SED suggests that it is 
low mass (^0.3 M 0 ), i.e., mid M. It powers a well-collimated 
bipolar molecular outflow (Mitchell et al. 1997). Compared to 
DG Tau the [O i] 63 /xm emission is clearly extended and much 
stronger in the red- shifted outflow lobe than on the star (Podio 
et al. 2012). At 70 /xm the continuum emission is unresolved 
on both stars, but at 100 /xm and especially 160 /xm one can 
see faint emission from the surrounding dark cloud, possibly 
from dust being heated and compressed by the outflows that are 
powered by DG Tau and DG Tau B. The PACS line spectroscopy 
from both stars is discussed in detail by Podio et al. (2012). 

DL Tau, a CTTS of spectral type K7 V, is classified as an 
outflow/jet source based only on the broad [O i] A6300 and [S n] 
A6371 (Hartigan et al. 1995). The [O i] A6300 line shows a high 
velocity component, which is blue-shifted to 300 km s _1 , 
while the [Sn] A6371 goes to —250 km s -1 , and has no low 
velocity component at all. DL Tau has deep and broad blue- 
shifted He i A 10830 (Edwards et al. 2003, 2006), which Kwan 
& Fischer (2011) show originates in an outflow rather than 
from accretion. C. Grady (2013, private communication) reports 
that a jet as well as two distant HH objects have recently been 
discovered in DL Tau. The disk is seen relatively face-on with 
an inclination i = 38 zb 2° (Guilloteau et al. 201 1). We find the 
[O i] A63 /xm emission to be in the outflow regime (Figure 6), 
which confirms that DL Tau indeed powers an outflow. 

GG Tau is a hierarchical quadruple system with a 10" sepa- 
ration between the two binaries (Leinert et al. 1991). GG Tau 
was one of the brightest millimeter sources in the survey by 
Beckwith et al. (1990) suggesting that it must be surrounded by 
a massive circumstellar disk. Millimeter interferometry (Simon 
& Guilloteau 1992; Guilloteau et al. 1999) shows that GG Tau 
Aa/Ab is surrounded by a large circumbinary disk with an outer 
radius of 2" -4". GG Tau Aa has a spectral type of K7, while 
the secondary, GG Tau Ab, has a spectral type of M0. 5 (White 
et al. 1999). The inner part of the ring/disk has been resolved in 
high-resolution scattered light images from 0.5 to 4 /xm (Rod- 
dier et al. 1996; McCabe et al. 2002; Duchene et al. 2004; Krist 



Right Ascension (2000.0) 

Figure 13. Logarithmically stretched 70 /xm image of GG Tau showing the 
trifoil structure of the PSF. The additional spur to the south coincides with 
the southern binary, GG Tau Bab, which is shown with star symbols. A PSF 
subtraction (not shown) confirms the detection of GG Tau B. 

(A color version of this figure is available in the online journal.) 

et al. 2005) and found to be elliptical with a semimajor axis 
of 1"5. 

The secondary pair, GG Tau Ba and Bb is also a binary and 
separated from the primary by ICE 4 at a PA. of 185° (Kraus & 
Hillenbrand 2009). GG Tau Ba has a spectral type of M5 and 
Bb M7 (White et al. 1999). We detect faint emission from the 
system at 70 /xm (Figure 13). 

HV Tau and DO Tau are two CTTSs associated with far- 
infrared nebulosities. The separation between the two is ~90'. / 8 
and both were imaged simultaneously with the PACS imager, 
although DO Tau is close to the edge of the image. Both stars are 
associated with far-infrared nebulosities, which become visible 
at 100 /xm (Figure 14). At 160 /xm the emission is very extended 
and connect both stars in a common envelope. 

DO Tau is a G star associated with an arc-like nebula, 
which is aligned with the nebulosity we see at 100 and 
160 /xm (McGroarty & Ray 2004). It drives a bipolar jet 
with the redshifted jet at a PA. of 70° (Hirth et al. 1994), 
approximately aligned with the arc-like nebulosity we see at 
100 /xm. McGroarty & Ray (2004) found three HH objects 
at P.A.s of 74°-78° northeast of DO Tau. The farthest one, 
HH 831 B, is IE from DO Tau. 

HV Tau is a close binary (Simon et al. 1996), which show 
no infrared excess (Woitas & Leinert 1998). About 4" to the 
northeast is a third star, HV Tau C, which has an edge on 
disk (Monin & Bouvier 2000; Terada et al. 2007). The spectral 
type of HV Tau C is K6 (White & Hillenbrand 2004). This 
is the star we see in the far-infrared, and which powers a 
bipolar jet (Stapelfeldt et al. 2003; McGroarty & Ray 2004; 
Duchene et al. 2010) with the PA. of the blue- shifted jet being 
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Right Ascension (2000.0) 

Figure 14. False color images of HV Tau and DO Tau at 70 /zm, 100 /zm, and 
160 jtzm overlaid with logarithmic contours. DO Tau is close to the edge of the 
image, which results in excess noise east of the star. Both stars appear to be 
associated with extended emission, which is barely visible at 70 /zm. At 160 /zm 
the emission is very extended. 

(A color version of this figure is available in the online journal.) 


~25° (McGroarty & Ray 2004). In this case the V-shaped 
100 /xm emission is aligned with a diffuse reflection nebula 
west-northwest of the star (Martin et al. 1994), which appears 
to be approximately aligned with the disk plane, and orthogonal 
to the jet. 

We find that both stars have strong [O i] 63 /xm emission, as 
expected from stars powering jets. 

HL Tau and XZ Tau are two of the first T Tauri stars found to 
power well-collimated optical jets (Mundt & Fried 1983; Mundt 
et al. 1990; Movsessian et al. 2007). HL Tau is a Class I object 
of spectral type K7 (White & Ghez 2001). It has very strong 
high velocity blue-shifted [O i] X 6300 emission (Edwards et al. 
1987) and it is one of the brightest T Tauri stars at 1.3 mm 
(Beckwith et al. 1990). The disk has an inclination angle of 
~42° (Lay et al. 1997). The [O i] 63 /xm emission is definitely in 


the outflow regime (Figure 6), but it does not have a particularly 
strong excess. XZ Tau, which has a lower accretion rate, has a 
much larger excess in [O i] 63 /xm due to the outflow. 

RY Tau is classified as a F8-G1 type star (Mora et al. 
2001; Calvet et al. 2004), although Kenyon & Hartmann (1995) 
assigned it a spectral type as late as K1 . It is a well known outflow 
source, which drives a well-collimated jet extending out at least 
31" from the star and with a counter-jet extending out to at least 
3'.5 in the opposite direction (St-Onge & Bastien 2008). The 
star has relatively strong free-free emission (Rodmann et al. 
2006), which is consistent with a thermal wind. The spectro- 
imaging observations by Agra-Amboage et al. (2009) detect the 
blue-shifted [O i] A6300 within 2" of the star, consistent with 
the large scale jet seen by St-Onge & Bastien (2008). Agra- 
Amboage et al. (2009) estimate the inclination of the jet to 
be between 45°-77°. Isella et al. (2010) resolved RY Tau at 
1.3 mm with millimeter interferometry and found that it has a 
large central cavity (28 AU) similar to what is seen in transitional 
disks. Alternatively the dust grains in the inner disk have grown 
to centimeter sizes lowering the dust opacity. Espaillat et al. 
(2011) found that they could fit the SED of RY Tau with an 
1 8 AU gap which contains some optically thin dust consistent 
with the cavity observed by Isella et al. (2010). 

RY Tau is also peculiar in other respects. Lommen et al. 
(2010) find a correlation between the strength of the 10 /xm sil- 
icate feature and the slope of the millimeter emission measured 
between 3 and 1 mm, except for RY Tau, which in their data 
set is an extreme outlier, although this is not entirely true. Our 
analysis, based on all published millimeter/sub-millimeter data, 
shows that the millimeter slope is rather normal, therefore in 
that respect RY Tau is not that unusual. The strong silicate emis- 
sion, however, requires an abundance of small grains. What is 
unique though is that a transitional disk star would drive a well- 
collimated optical jet, have free-free emission and still be rather 
strongly accreting (M ~ 6-9 x 10 -8 M 0 yr -1 ; Agra-Amboage 
et al. 2009). 

Our observations shows no excess from the jet in the [O i] 
63 /xm line. The line intensity is completely consistent with a 
disk origin (Figure 6). This could mean either that there is no 
[O i] 63 /xm emission from the outflow or that the emission from 
the disk is anomalously low, perhaps due to the central cavity. 

StHA 34 (HBC 425) is a spectroscopic binary with compo- 
nents of nearly equal luminosity and temperature (both M3) and 
broad, strong Ha emission characteristic of a CTTS (White & 
Hillenbrand 2005). However, neither component of the binary 
shows Li i A 6708 absorption suggesting a very long-lived accre- 
tion disk. Comparison with PMS evolutionary models give an 
isochronal age of 8 zb 3 Myr, which is much younger than the 
predicted lithium-depletion timescale of ~25 Myr. Hartmann 
et al. (2005), who observed StHA 34 with IRS on Spitzer mod- 
eled the infrared SED with three components: an inner disk wall, 
an optically thin inner disk, and an outer disk, i.e., similar to 
the SED of a transitional disk. They favored an age of 25 Myr, 
which would require StHA34 to be at a distance of ~100 pc, 
and therefore not a member of the Taurus association. 

Dahm & Lyke (2011) found a low-mass companion ~ 1^23 
southeast of the primary pair, which has strong Lii A6708 
absorption. STHA 34 C has a spectral type of M5.5. Comparison 
with PMS evolutionary tracks imply a mass of ~0.09 M© and 
an age of 8-10 Myr assuming the nominal distance to Taurus 
(140 pc). It therefore seems more likely that StHA 34 is lithium 
depleted and not surrounded with a 25 Myr old, still accreting 
circumbinary disk. 
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We detected StHA 34 in continuum at 70, 100, and 160 /x m 
(Table 1). The emission is rather faint, which is consistent with 
an evolved, possibly transitional disk. 

T Tauri is an exceptional member of the class of PMS stars 
that share its name. It is a hierarchical triple system with the 
optical star T Tau N having a K-type spectrum and large optical 
variations on timescales of months to years (Beck & Simon 
2001), and the southern deeply embedded binary system T Tau 
S, which has a projected separation of 0?05 (Koresko 2000). 
T Tau N is known to power a bright optical jet (Biihrke et al. 
1986). T Tau S is surrounded by a massive accretion disk and 
powers a spectacular jet (Reipurth et al. 1997). 

ISO found T Tauri to have the richest line mid/FIR emission 
line spectrum of any PMS star (Lorenzetti 2005). T Tauri is very 
bright in the far-infrared (Table 2), and the emission is extended 
(~4") in all PACS bands. We see faint emission from the Hind’s 
nebula at 70 /xm, which is much stronger at 100 and 160 /x m. At 
160 /xm we also start to pick up emission from the surrounding 
cloud. 

Since T Tauri is a triple system which drives at least two 
optical jets, it is not clear where the emission originates. It 
could arise in the disk, in the outflow, or in the hot surrounding 
envelope (Podio et al. 2012). We have therefore excluded it from 
any correlation analysis. 

V773 Tau is a compact quadruple system with at least 
four stars within 0'. / 3 (Boden et al. 2012). V773 Tau A is a 
spectroscopic binary with spectral types of K2 and K5 (Welty 
1995). It is a WTTS and one of the strongest radio stars in Taurus 
(O’Neal et al. 1990; Massi et al. 2008). The radio emission is 
non-thermal and strongly variable (Dutrey et al. 1996). The 
variability is due to interacting coronae, which causes strong 
flaring at periastron (Massi et al. 2008). V773 Tau A has no 
near infrared excess, and is therefore unlikely to be surrounded 
by a circumbinary disk. V773 Tau B is a CTTS, which varies 
by more than 0.5 mag in the visible and has a spectral type of 
K7-M0.5 (White & Ghez 2001; Duchene et al. 2003; Boden 
et al. 2012). It has some near-IR excess longward of 2 /x m and 
maybe surrounded by a compact disk. V773 Tau C is extremely 
red and almost certainly surrounded by a circumstellar disk 
(Duchene et al. 2003). Since V773 Tau was detected at 1.3 mm 
and 850 /xm (Osterloh & Beckwith 1995; Andrews & Williams 
2005) it is clear that the system must have a circumstellar 
disk. Furthermore we easily detected it at both 60 and 160 /xm 
(Table 2). Whether the dust emission is dominated by V773 Tau 
C or whether V773 Tau B also contributes requires observations 
with sub-arcsecond imaging resolution with good sensitivity, as 
can be provided by ALMA. 

V773 Tau is not known to drive a jet, yet the [O i] 63 /x m flux 
density clearly puts it in the outflow regime (Figure 6). Cabrit 
et al. (1990) did detect [O i] A6300 and [S n] A673 1 emission but 
did not publish a spectrum. It is not clear whether the emission 
comes from a jet, especially since the velocity offsets were small. 
However, FM Tau, an M0 star, is 38" to the north-northwest 
of V773 Tau, and CW Tau, a K3 star, 92" to southeast. The 
latter is a young Class II type CTTS, which powers a bipolar 
optical jet and excites several HH objects (Dougados et al. 2000; 
McGroarty & Ray 2004). Not all the HH objects discovered by 
McGroarty & Ray (2004) can be explained by the outflow from 
CW Tau. In particular, HH 828, which is just north of FM Tau 
and V773 Tau, is not aligned with the CW Tau jet, and therefore 
more likely to be powered either by FM Tau or V773 Tau. The 
[O i] X 6300 line in FM Tau is slightly blue-shifted and does not 
show much high velocity emission (Hartigan et al. 1995), and 


is therefore an unlikely source. V773 Tau may well be exciting 
HH 828. In our analysis we therefore assume that V773 Tau is 
an outflow source, although this still needs to be confirmed. 

6. SUMMARY AND CONCLUSIONS 

In this paper we present the results of Herschel / PACS 
photometric and spectroscopic observations of a large sample 
of T Tauri stars in the Taurus/ Auriga star-forming regions as 
part of the GASPS survey. We observed 82 fields in photometry 
and 75 in spectroscopy, for a total sample of about 120 objects, 
mostly Class II/III objects (CTTS /WTTS) as well as a handful 
of embedded Class I sources and transitional disks. 

In photometry, we detected 50 known members of the 
association, apart for 3 Class III sources (FW Tau, StHA 34, 
and V819 Tau), all of them are Class II or Class I sources. Most 
objects are spatially unresolved in the continuum. Exceptions 
include T Tau, the multiple systems UZ Tau and GG Tau and 
HD 283759 (but the latter is not a member of the association). 

Even though our spectroscopic observations include a slew of 
other molecular transitions, we only discuss here the results for 
[O i] 63 /xm and 145 /x m, as well as [C n] 158 /xm. Results on the 
0 -H 2 O line at 63.32 /xm were presented in Riviere-Marichalar 
et al. (2012) while other transitions, typically detected in 
strong outflow sources and sometimes spatially resolved, were 
discussed in Podio et al. (2012) for a subset of our sample. 

The [0 1 ] 63 /xm line, spatially and spectrally unresolved in 
most sources, is by far the strongest line in our sample and often 
the only line detected, especially for fainter sources. With a 
typical line sensitivity of ~1 x 10“ 17 W m -2 , we obtained >3 a 
detections in 43 out of 75 targets; no Class III sources were 
detected. Marginal (2a -3a level) detections are presented for 
another 10 sources. The [0 1 ] 145 /xm line is 10-25 times fainter 
than the [Oi] 63 /xm line and detected only in 17 sources. The 
detection rate for [Cn] 158 /xm line is even lower (10 targets) 
and detections are restricted to jet/outflow sources. 

We found no correlation between [0 1 ] 63 /xm line intensity 
and disk mass, suggesting that the line is either optically thick or 
is not a valuable direct probe the mass of the overall gas reservoir. 
However, we found a tight correlation between the [0 1 ] 63 /xm 
line strength and the 63 /xm continuum flux density for “disk- 
only” objects, i.e., sources with no or only weak outflow activity. 
Sources with clear evidence for outflow activity, on the other 
hand, can show line strengths that are as much as 20 times 
higher than disk-only sources, suggesting that the line emission 
is dominated by the jet. Some outflow sources (6 out of 24), 
however, follow the disk correlation and show little or no excess 
due to the outflow. Among disk sources, transitional disks also 
stand out, being fainter in [0 1 ] than CTTSs by a factor of two 
to three for comparable continuum flux densities. 

The tight correlation between [0 1 ] line and 63 /xm continuum 
emissions in the disk-only sample suggests that they originate 
from the same region of the disk. Simple disk models indicate 
that the most of the far-infrared continuum emission arise from 
the inner disk regions (within 5-50 AU at most, depending 
on the disk outer radius). Based on two-component isothermal 
graybody fits to the millimeter and far-infrared SEDs, we 
suggest that the [0 1 ] 63 /xm emission and the 63 /xm continuum 
emission originate from the inner disk and/or lower surface 
layers where the gas and the dust is thermalized. This can 
also explain the weaker line emission of transition disks, since 
the dust in these systems is substantially cooler. The absence 
of dependence of both the [0 1 ] 63 /xm line and far-infrared 
continuum emission strength on binary separation for systems 
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wider than 10 AU further supports the co-spatiality of their 
emitting regions and the fact that they are limited to the 
innermost regions of the disk. 

The observations presented in this paper provide a rich dataset 
to test thermo-chemical models which will provide a better 
understanding of the origin and strength of the [O i] 63 /xm 
line in protoplanetary disks. 
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